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Abstract 


NIST  has  completed  Phase  n.A  of  a project  to  study  the  impact  of  HVAC  systems  on  residential 
indoor  air  quality  and  to  assess  the  potential  for  using  residential  forced-air  systems  to  control 
indoor  pollutant  levels.  In  this  effort,  NIST  is  performing  whole  building  airflow  and 
contaminant  dispersal  computer  simulations  with  the  program  CONTAM93  to  assess  the  ability 
of  modifications  of  central  forced-air  heating  and  cooling  systems  to  control  pollutant  sources 
relevant  to  the  residential  environment.  This  report  summarizes  the  results  of  Phase  n.A  of  this 
project,  which  consisted  of  three  major  efforts:  baseline  simulations  of  contaminant  levels 
without  indoor  air  quality  (lAQ)  controls,  design  of  the  lAQ  control  retrofits,  and  preliminary 
simulations  of  contaminant  levels  with  the  lAQ  control  retrofits.  In  Phase  H.B  of  the  study,  all  of 
the  baseline  cases  will  be  modified  to  incorporate  the  lAQ  control  retrofits.  The  retrofit  results 
will  then  be  compared  to  the  baseline  results  to  evaluate  the  effectiveness  of  the  retrofits. 

The  pollutant  concentrations  in  a building  depend  on  many  factors  including  the  configuration  of 
the  building  zones,  the  air  leakage  of  the  building  envelope  and  of  interior  partitions,  wind 
pressure  profile  on  the  building  envelope,  pollutant  source  strengths  and  temporal  profiles, 
heating  and  cooling  system  airflow  rates,  furnace  filter  efficiency,  characteristics  of  reversible 
pollutant  sinks  in  the  building,  individual  pollutant  decay  or  deposition  rates,  and  ambient 
weather  and  pollutant  concentrations.  This  report  describes  the  input  data  used  to  model  the 
baseline  houses  with  CONTAM93  and  presents  the  results  of  the  baseline  simulations  in  the 
form  of  the  transient  pollutant  concentrations  for  selected  simulations  and  a summary  of  peak 
and  average  concentrations  for  all  baseline  simulations.  Three  indoor  air  quality  control 
technologies  were  then  selected  for  incorporation  into  the  baseline  house  models  to  determine 
their  effectiveness  in  controlling  the  modeled  pollutant  sources.  The  technologies  include  the 
following:  electrostatic  particulate  filtration,  heat  recovery  ventilation,  and  an  outdoor  air  intake 
damper  on  the  forced-air  system  return.  Selected  baseline  cases  were  then  modified  to  implement 
these  indoor  air  quality  control  retrofits,  and  preliminary  simulations  were  performed  to 
demonstrate  the  ability  of  the  program  to  model  the  control  technologies. 


Key  Words:  airflow  modeling,  building  technology,  computer  simulation,  filtration,  heat 
recovery  ventilator,  HVAC  system,  indoor  air  quality,  infiltration,  residential  buildings, 
ventilation 
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Introduction 


Despite  the  increasing  interest  in  residential  indoor  air  quality  (lAQ)  problems,  only  limited 
research  has  been  conducted  which  integrates  the  analysis  of  pollutant  sources,  residential 
heating  and  cooling  system  operation,  and  building  characteristics.  While  central  forced-air 
heating  and  cooling  systems  may  provide  solutions  to  some  lAQ  problems,  such  an  integrated 
approach  is  required  to  analyze  these  options.  Because  large  quantities  of  indoor  air  circulate 
through  these  systems,  they  offer  the  potential  for  treating  the  indoor  air  and  then  distributing  this 
treated  air  through  the  system  ductwork  to  the  building.  Also,  outdoor  air  brought  into  the 
building  by  the  forced-air  system  can  be  distributed  throughout  the  building  by  this  ductwork. 
Most  modeling  studies  of  lAQ  in  residential  buildings  have  employed  very  simple  models  of  the 
building  and  its  systems,  ignoring  the  multizone  nature  of  the  airflows  involved.  The  use  of  such 
simple  analytical  procedures  has  limited  our  understanding  of  the  impact  of  central  forced-air 
heating  and  cooling  systems  on  residential  lAQ  and  the  possibility  of  using  these  and  other 
systems  to  mitigate  lAQ  problems. 

The  National  Institute  of  Standards  and  Technology  (NIST)  is  conducting  a study  for  the  U.S. 
Consumer  Product  Safety  Commission  (CPSC)  to  assess  the  potential  effectiveness  of  existing 
heating,  ventilating,  and  air  conditioning  (HVAC)  technology  to  reduce  the  levels  of  selected 
pollutants  in  single-family  residential  buildings.  This  effort  is  employing  a new  multizone 
airflow  and  contaminant  dispersal  modeling  program,  CONTAM93  (1).  In  this  effort,  NIST  is 
performing  whole  building  airflow  and  contaminant  dispersal  computer  simulations  to  assess  the 
ability  of  modifications  of  central  forced-air  heating  and  cooling  systems  to  control  pollutant 
sources  relevant  to  the  residential  environment.  Phase  I of  the  project  included  conducting  a 
literature  review,  developing  a detailed  simulation  plan,  and  hosting  a workshop  to  discuss  the 
project,  and  was  described  in  a previous  report  (2).  This  report  summarizes  the  results  of  Phase 
11.  A which  consisted  of  three  major  tasks:  baseline  simulations  of  contaminant  levels  without 
lAQ  controls,  design  of  the  LAQ  control  retrofits,  and  preliminary  simulations  of  contaminant 
levels  with  the  lAQ  control  retrofits  in  place. 

In  Phase  n.B  of  the  study,  the  baseline  HVAC  systems  will  be  modified  to  incorporate  the  lAQ 
control  technologies  described  in  this  report  and  simulations  will  be  performed  for  all  conditions 
under  which  baseline  simulations  were  performed.  The  Phase  II.B  simulation  results  will  be 
compared  with  the  results  presented  here  to  determine  the  effectiveness  of  the  lAQ  modifications 
at  controlling  the  selected  pollutant  sources. 

Contents  of  Report 

The  first  section  of  the  report  describes  the  baseline  simulations  performed.  The  program 
CONTAM93  (1)  was  used  to  calculate  airflows  and  pollutant  distributions  for  the  houses  and 
pollutant  sources  described  in  the  report  on  Phase  I of  the  project  (2).  The  houses  modeled  are 
not  based  on  real  buildings  but  are  intended  to  be  representative  of  typical  buildings.  This  first 
section  presents  the  input  data  used  to  describe  the  houses,  HVAC  systems,  pollutants,  sources, 
and  boundary  conditions  in  the  baseline  simulations.  In  addition,  this  section  summarizes  the 
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results  of  the  baseline  simulations  including  transient  pollutant  concentrations  for  selected 
simulations  and  a summary  of  peak  and  average  concentrations  for  all  baseline  simulations. 

The  second  section  describes  the  indoor  air  quality  control  technologies  that  will  be  evaluated  in 
the  computer  simulations  during  Phase  II.B.  These  technologies  will  be  incorporated  into  the 
baseline  house  models  to  determine  their  effectiveness  in  controlling  the  selected  pollutant 
sources.  The  three  technologies  described  in  this  section  include  electrostatic  particulate 
filtration,  heat  recovery  ventilation,  and  an  outdoor  air  intake  damper  on  the  forced-air  system 
return.  This  section  describes  each  of  these  technologies  and  includes  revisions  of  the  baseline 
house  duct  drawings.  In  addition,  this  section  contains  an  estimate  of  the  equipment  and 
installation  costs  and  a revision  of  the  thermal  load  calculations  based  on  the  modifications. 
Finally,  the  impacts  of  each  of  these  technologies  on  “other  contaminants”  are  discussed 
qualitatively.  These  other  contaminants,  as  described  in  the  original  project  work  statement, 
include  contaminants  that  have  typically  been  of  concern  to  designers  of  residential  ventilation 
systems  including  cooking  odors,  tobacco  smoke,  moisture,  outdoor  pollen,  outdoor  odors  and 
ozone. 

The  third  section  presents  the  results  of  preliminary  simulations  of  the  lAQ  control  retrofits. 
These  simulations  involved  modifying  selected  baseline  simulation  cases  with  the  three  lAQ 
control  retrofits.  The  preliminary  simulations  were  performed  to  demonstrate  the  ability  of  the 
program  to  model  the  lAQ  control  technologies. 

The  report  includes  two  appendices.  The  first  appendix  describes  modeling  performed  to 
characterize  the  airflow  in  the  houses  including  the  results  of  fan  pressurization  simulations  and 
whole  house  infiltration  simulations.  The  second  appendix  includes  summary  tables  of  the 
baseline  and  preliminary  simulation  results. 
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Baseline  Simulations 


This  section  of  the  report  describes  the  baseline  simulations  performed  in  Phase  11.  A.  This 
section  presents  the  input  data  describing  the  houses,  HVAC  systems,  pollutants,  sources,  and 
boundary  conditions  modeled  in  the  baseline  simulations.  In  addition,  this  section  summarizes 
the  results  of  the  baseline  simulations  including  transient  pollutant  concentrations  for  selected 
simulations  and  a summary  of  peak  and  average  concentrations  for  all  baseline  simulations. 

Baseline  Simulation  Input  Data 

Calculating  airflow  rates  and  contaminant  concentrations  with  CONTAM93  or  any  other 
multizone  model  requires  the  following  input:  the  configuration  and  volume  of  the  building 
zones,  the  air  leakage  paths  through  the  building  envelope  and  interior  walls,  wind  pressure 
profile  on  the  building  envelope,  pollutant  source  strengths  and  temporal  profiles,  HVAC  system 
flows,  furnace  filter  efficiency,  characteristics  of  reversible  pollutant  sinks,  individual  pollutant 
decay  or  deposition  rates,  and  ambient  weather  and  pollutant  concentrations.  This  section 
describes  the  input  data  used  in  the  baseline  simulations. 

The  study  included  eight  building  models  - a ranch  and  a two-story  house,  located  in  two  sites 
(Miami  and  Minneapolis),  with  typical  and  low  values  of  airtightness.  The  Phase  I NISTIR  (2) 
described  the  layout  and  dimensions  of  each  house  and  contained  floorplan  drawings. 
Simulations  were  performed  under  three  sets  of  weather  conditions  (cold,  mild,  and  hot)  for  each 
building.  Each  simulation  was  performed  for  a one-day  cycle  repeated  until  peak  concentrations 
converged  to  a specified  tolerance.  Referring  to  all  pollutant  sources  modeled  for  a single 
building  as  one  simulation,  there  were  a total  of  24  baseline  simulation  cases.  Table  1 lists  the 
baseline  simulations  by  house  type,  location,  airtightness  and  weather  condition. 
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Table  1 - Baseline  simulations 


Simulation 

House  type 

Location 

Airtightness 

Weather 

SIMIFLC 

ranch 

Miami 

typical 

cold 

SIMIFLM 

ranch 

Miami 

typical 

mild 

SIMIFLH 

ranch 

Miami 

typical 

hot 

SIMIFTC 

ranch 

Miami 

tight 

cold 

SIMIFTM 

ranch 

Miami 

tight 

mild 

SIMIFTH 

ranch 

Miami 

tight 

hot 

SIMIMLC 

ranch 

Minneapolis 

typical 

cold 

SIMIMLM 

ranch 

Miimeapolis 

typical 

mild 

SIMIMLH 

ranch 

Minneapolis 

typical 

hot 

SIMIMTC 

ranch 

Minneapolis 

tight 

cold 

SIMIMTM 

ranch 

Minneapolis 

tight 

mild 

SIMIMTH 

ranch 

Minneapolis 

tight 

hot 

SIM2FLC 

two-story 

Miami 

typical 

cold 

SIM2FLM 

two-story 

Miami 

typical 

mild 

SIM2FLH 

two-story 

Miami 

typical 

hot 

SIM2FTC 

two-story 

Miami 

tight 

cold 

SIM2FTM 

two-story 

Miami 

tight 

mild 

SIM2FTH 

two-story 

Miami 

tight 

hot 

SIM2MLC 

two-story 

Minneapolis 

typical 

cold 

SIM2MLM 

two-story 

Minneapolis 

typical 

mild 

SIM2MLH 

two-story 

Minneapolis 

typical 

hot 

SIM2MTC 

two-story 

Minneapolis 

tight 

cold 

SIM2MTM 

two-story 

Minneapolis 

tight 

mild 

SIM2MTH 

two-story 

Minneapolis 

tight 

hot 

Detailed  information  on  building  component  leakage  of  the  houses  is  not  available  as  the  houses 
modeled  were  not  based  on  real  buildings.  However,  since  there  is  no  attempt  to  compare 
predictions  with  experimental  data,  the  building  leakage  modeled  needs  only  to  be  reasonable  in 
magnitude  and  distribution.  Table  2 shows  all  of  the  leakage  paths  between  the  zones  of  the 
Miami  ranch  house  (see  Figure  1 for  the  ranch  house  floorplan  and  zone  labeling  and  Figure  2 for 
the  two-story  house  floorplan  and  zone  labeling).  The  Minneapolis  houses  have  basements  (zone 
label  BMT)  that  are  not  shown  in  the  figures.  Table  3 lists  the  values  for  those  leakage  paths  for 
both  the  typical  and  tight  cases.  The  Table  3 leakage  areas  are  for  a reference  pressure  difference 
of  4 Pa  and  a discharge  coefficient  of  1.0  and  are  based  on  values  listed  in  Table  23-3  of 
ASHRAE  (3)  unless  otherwise  noted.  The  typical  values  were  generally  based  on  "best  estimate" 
and/or  imcaulked  entries  in  the  ASHRAE  table,  while  the  tight  values  were  based  on  minimum 
and/or  caulked  entries.  All  doors  connecting  zones  other  than  closets  were  modeled  as  open.  The 
same  leakage  values  were  used  for  the  other  houses,  although  the  paths  connecting  the  zones 
differed  depending  on  the  house  configurations. 
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Figure  1 - Ranch  house  floorplan  and  zones 


Figure  2 - Two-story  house  floorplan  and  zones 
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INTW 

OUTL 


INTW 

OUTL 


MBA^ 


INTD 

INTW 
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Table  2 - Air  leakage  paths  for  Miami  ranch  house 


INTW 

OUTL 


INTW 

OuTi. 


INTW 


INTW 


INTW 

OUTL 


INTW 

OUTL 


INTW 

OUTL 


INTW 

INTD 

OUTL 


INTD 

INTW 


INTD 

INTW 


INTD 

INTW 

OUTL 


INTD 

INTW 


OLD 

INTW 


INTD 

INTW 

OUTL 


HAD 


EXTD 

EXW 

OUTL 


CEIL 

CPEN 


CEIL 

CPEN 


CEIL 

CPEN 


CEIL 

CPEN 

PIP 


CEIL 

CPEN 

PIP 


CEIL 

CPEN 


CEIL 

CPEN 


CEIL 

CPEN 

CPEN 


CEIL 

ATD 


EXV 

EXW 

OUTL 


WIN 

EXW 

OUTL 


WIN 

EXW 

OUTL 


WIN 

EXW 

OUTL 


EXV 


WIN 

EXV 

EXW 

OUTL 


SGD 

EXTD 

WIN 

EXW 

OUTL 


GAD 

GARF 

EXW 


VNT 


The  infiltration  through  a building's  envelope  also  depends  on  the  static  pressure  distribution 
created  by  the  wind  on  the  building's  exterior  surfaces.  The  relationship  between  wind  and 
surface  pressures  are  characterized  by  wind  pressure  coefficients  which  depend  on  the  wind 
direction,  the  building  shape,  the  position  on  the  building  surface,  and  the  presence  of  shielding 
near  the  building.  The  surface  pressure  coefficients  for  the  building  walls  were  based  on 
Equation  23-8  of  ASHRAE  (3).  The  coefficient  for  the  flat  garage  roof  was  based  on  Figure  14-6 
of  ASHRAE  (3).  The  ASHRAE  wind  pressure  coefficients  do  not  include  shielding  effects  and 
no  modifier  for  shielding  effects  was  used,  however,  recent  studies  have  reported  on  the 
shielding  effects  of  trees  (4)  and  rows  of  houses  (5). 
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Table  3 - Air  leakage  values 


Name 

Description 

Typical 

Tight 

ATD 

Attic  door 

30  cm^/ea 

18  cmVea 

CEIL 

Ceiling  [Based  on  general  ceiling] 

1.8  cmVm^ 

0.79  cmVm' 

CLD 

Closet  door  (closed)  [Based  on  interior  door] 

0.9  cm^/m 

0.25  cmVm 

Closet  door  frame  [Based  on  general  door  frame] 

25  cmVea 

12  cmVea 

CPEN 

HVAC  ceiling  penetration  [Based  on  kitchen  vent  with  damper  closed] 

5 cmVea 

1 cmVea 

EXTD 

Exterior  door  [Single] 

21  cmVea 

12  cmVea 

Door  frame  [Wood] 

1.7  cmVm^ 

0.3  cmVm' 

EXV 

Bathroom  exhaust  vent 

20  cmVea 

10  cmVea 

Kitchen  exhaust  vent 

40  cmVea 

5 cm'/ea 

EXW 

Ceiling- wall  joint 

1 .5  m^/m 

0.5  mVm 

Floor- wall  joint 

4 cm^/m 

0.8  cmVm 

Wall-wall  joint  [Based  on  ceiling-wall  joint] 

1.5  mVm 

0.5  mVm 

GAD 

Garage  door  [Based  on  general  door  (2  m x 4 m)] 

0.45  cmVm 

0.31  cmVm 

Garage  door  frame  [Wood] 

1.7  cmVm^ 

0.3  cmVm' 

GARF 

Garage  roof  [Based  on  general  ceiling] 

1.8  cmVm^ 

0.79  cmVm' 

HAD 

Hall  doorway 

2.4  mVea 

2.4  mVea 

INTD 

Interior  door  (closed)  [Based  on  Table  4.2  of  Klote  and  Milke  (6)] 

140  cmVea 

75  cmVea 

Interior  door  (open) 

2.1  mVea 

2.1  mVea 

INTW 

Interior  wall  [Based  on  gypsum  board  on  stud  wall  (Shaw  et  al.  7)] 

2.0  cm^/m^ 

2.0  cmVm' 

OUTL 

Electric  outlet 

2.5  cmVea 

0.5  cmVea 

PIP 

Piping  penetrations 

6 cm^/ea 

2 cmVea 

SGD 

Sliding  glass  door 

22  cmVea 

3 cm'/ea 

VNT 

Attic  vent  [Based  on  Table  21-1  of  3] 

1 cmV300  cm' 

1 cm'  / 300  cm' 

WIN 

Double  hung  window 

2.5  cmVm 

0.65  cm'/m 

Window  framing  [Wood] 

1.7  cmVm' 

0.3  cm'/m' 

The  building  HVAC  systems  were  designed  in  Phase  I of  the  study  and  are  described  in  the  Phase 
I report  (2).  This  earlier  report  contains  the  heating  and  cooling  equipment  types  and 
descriptions,  overall  and  individual  supply  and  return  airflow  rate  design  values  for  both  heating 
and  cooling,  and  drawings  showing  the  system  equipment  and  duct  locations  and  duct  sizes.  In 
addition  to  this  information,  CONTAM93  requires  information  on  the  system  operation, 
specifically,  an  on-off  schedule.  The  schedule  was  determined  by  calculating  the  fractional 
on-time  required  to  meet  the  cooling  or  heating  load  for  each  3 -hour  period  of  the  day.  A control 
profile  incorporating  this  schedule  was  then  input  for  each  simulation.  For  the  baseline 
simulations,  the  HVAC  systems  included  standard  furnace  filters  with  constant  efficiencies  of 
5%  for  fine  particles  and  90%  for  coarse  particles.  Fine  particles  are  defined  as  having  a diameter 
less  than  2.5  fxm  (the  efficiency  is  based  on  a diameter  of  0.6  jim);  coarse  particles  are  defined  as 
having  a diameter  greater  than  2.5  |im  (the  efficiency  is  based  on  a diameter  of  6 |xm).  These 
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efficiency  values  are  based  on  assumed  arrestance  for  these  filters  of  about  90%  and  a review  of 
manufacturers'  test  data.  No  outdoor  air  intake  is  included  for  the  baseline  HVAC  systems. 

Another  important  consideration  for  the  HVAC  systems  is  duct  leakage.  Since  the  duct  system 
itself  is  not  modeled  in  these  simulations,  duct  leakage  is  modeled  by  including  an  additional 
system  supply  or  return  point  and  reducing  the  other  supply  and  return  flows  by  the 
corresponding  amount  rather  than  by  an  effective  leakage  area.  Cummings  et  al.  1991  (8)  tested 
duct  leakage  in  160  houses  in  Florida  and  found  that  return  leaks  were  dominant  in  the  majority 
of  homes.  They  reported  an  average  return  leak  fraction  of  10.7%  (based  on  ratio  of  leakage  flow 
to  total  system  flow).  For  the  Minneapolis  houses,  a return  leak  of  10%  was  included  in  the 
basement.  For  the  Miami  ranch  house,  a supply  leak  of  10%  was  included  in  the  attic  because  the 
system  has  a central  return.  For  the  Miami  two-story  house,  no  leaks  were  included  because  all 
ducts  and  equipment  are  internal . The  air  distribution  system  layouts  were  designed  based  on 
guidelines  published  by  the  National  Association  of  Home  Builders  (9)  and  drawings  are 
included  in  the  Phase  I report  (2). 

The  ambient  boundary  conditions  required  by  CONTAM93  include  weather  and  outdoor 
pollutant  concentrations.  The  weather  conditions  were  chosen  by  selecting  a hot,  mild,  and  cold 
day  for  each  location  from  Weather  Year  for  Energy  Calculation  (WYEC)  data  (10).  The  WYEC 
data  is  presented  in  Tables  5 and  6 for  Miami  and  Minneapolis,  respectively,  and  includes 
temperature,  wind  speed,  and  wind  direction  from  north. 
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Table  5 - Miami  weather  data 


Cold 

Mild 

Hot 

Hour 

T 

V 

wind 

Dir 

T 

V 

wind 

Dir 

T 

V 

wind 

Dir 

CO 

(m/s) 

C) 

CO 

(m/s) 

C) 

CO 

(m/s) 

C) 

0 

2.8 

2.3 

320 

13.3 

3.9 

360 

26.7 

0 

0 

1 

2.8 

2.3 

300 

13.3 

2.7 

360 

26.1 

1.2 

200 

2 

2.8 

3.5 

310 

13.3 

3.5 

360 

26.1 

1.2 

200 

3 

2.8 

2.7 

320 

13.9 

1.9 

20 

25.6 

1.6 

200 

4 

2.2 

2.3 

310 

13.3 

2.7 

20 

25.6 

1.9 

200 

5 

2.2 

3.5 

310 

13.9 

1.6 

360 

26.1 

1.9 

230 

6 

2.8 

2.7 

320 

13.3 

2.3 

340 

25.6 

1.9 

200 

7 

3.3 

3.5 

300 

14.4 

2.3 

340 

26.7 

1.6 

230 

8 

4.4 

2.3 

290 

16.1 

2.7 

340 

27.2 

1.9 

200 

9 

6.1 

2.7 

330 

21.1 

4.7 

70 

30.6 

2.3 

200 

10 

8.9 

3.1 

320 

23.3 

4.7 

70 

31.7 

2.3 

230 

11 

11.7 

2.3 

320 

23.3 

5.1 

70 

32.8 

0 

0 

12 

13.9 

2.7 

330 

23.3 

5.4 

70 

33.3 

2.3 

200 

13 

14.4 

2.7 

350 

22.8 

5.1 

70 

33.3 

3.9 

140 

14 

16.1 

2.3 

360 

22.8 

5.4 

70 

32.8 

4.3 

180 

15 

17.2 

0.8 

40 

22.2 

4.7 

70 

31.7 

4.7 

160 

16 

17.8 

2.7 

40 

21.7 

3.9 

90 

30.6 

1.9 

290 

17 

17.2 

3.5 

20 

21.7 

3.1 

90 

31.7 

3.1 

140 

18 

16.7 

1.9 

340 

21.7 

4.3 

70 

30.6 

2.3 

160 

19 

16.1 

2.3 

340 

21.1 

4.3 

90 

27.8 

1.6 

50 

20 

15 

1.6 

350 

21.1 

2.7 

90 

27.8 

1.2 

50 

21 

14.4 

1.9 

350 

21.1 

3.1 

90 

27.2 

1.6 

200 

22 

16.1 

2.3 

30 

21.7 

1.2 

90 

26.1 

2.3 

230 

23 

16.1 

2.3 

60 

21.7 

2.3 

90 

26.1 

1.2 

250 

24 

17.2 

3.5 

60 

20.6 

3.1 

50 

26.1 

0 

0 
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Table  6 - Minneapolis  weather  data 


Cold 

Mild 

Hot 

Hour 

T 

(°C) 

V 

wind 

(m/s) 

Dir 

n 

T 

(°C) 

V 

wind 

(m/s) 

Dir 

n 

T 

CC) 

V 

wind 

(m/s) 

Dir 

n 

0 

-21.1 

1.6 

330 

7.8 

1.9 

60 

21.1 

3.1 

180 

1 

-21.1 

1.6 

330 

7.8 

1.9 

40 

20 

2.7 

180 

2 

-21.1 

3.1 

350 

7.8 

3.1 

90 

18.9 

2.7 

180 

3 

-21.1 

3.1 

350 

7.2 

1.9 

100 

17.8 

1.9 

180 

4 

-21.1 

3.1 

350 

7.2 

4.7 

130 

18.3 

1.6 

158 

5 

-21.1 

3.1 

350 

7.2 

3.9 

130 

17.2 

2.7 

135 

6 

-21.7 

3.5 

350 

7.2 

3.1 

120 

17.8 

3.5 

158 

7 

-21.7 

2.7 

340 

7.2 

3.9 

140 

20 

1.9 

158 

8 

-21.7 

2.7 

350 

7.8 

2.7 

120 

24.4 

4.7 

180 

9 

-21.1 

3.9 

340 

8.9 

3.1 

130 

26.1 

5.8 

180 

10 

-20.6 

3.9 

310 

7.8 

4.3 

130 

28.3 

6.6 

203 

11 

-20.6 

4.7 

310 

8.3 

4.7 

130 

30 

6.2 

203 

12 

-20.6 

3.9 

320 

8.9 

4.3 

140 

30.6 

6.2 

203 

13 

-20.6 

4.3 

320 

8.9 

4.7 

140 

31.1 

7 

203 

14 

-20 

5.1 

300 

8.3 

6.2 

120 

31.1 

7.4 

203 

15 

-20 

4.7 

290 

8.9 

6.2 

110 

31.1 

6.6 

203 

16 

-20.6 

4.3 

310 

8.9 

5.8 

130 

31.1 

6.6 

203 

17 

-21.1 

3.5 

290 

9.4 

5.1 

130 

28.9 

4.7 

203 

18 

-22.8 

3.1 

280 

9.4 

5.4 

130 

29.4 

4.7 

180 

19 

-23.3 

2.7 

280 

11.1 

5.4 

160 

27.8 

4.7 

180 

20 

-24.4 

3.1 

300 

11.7 

5.8 

170 

26.1 

4.3 

180 

21 

-25 

3.1 

280 

11.1 

6.2 

180 

24.4 

3.9 

180 

22 

-25.6 

2.7 

280 

11.1 

5.8 

200 

23.9 

3.9 

180 

23 

-27.2 

2.3 

240 

10.6 

6.2 

220 

23.3 

4.7 

158 

24 

-28.9 

2.3 

240 

7.8 

2.7 

240 

22.8 

4.3 

180 

Outdoor  pollutant  concentrations  vary  by  location  and  over  time  at  any  one  location.  The 
concentrations  used  as  boundary  conditions  for  the  indoor  sources  in  the  simulations  were 
selected  as  typical  outdoor  conditions  and  are  not  meant  to  represent  the  actual  conditions  at  any 
specific  location.  The  values  used  were  specified  per  the  schedules  in  Table  7.  The  CO  and  NO2 
concentrations  were  chosen  based  on  review  of  US  EPA  air  quality  documents  (11,  12,  13).  They 
were  chosen  to  have  a diurnal  pattern  with  morning  and  afternoon  peaks.  The  selected  CO  and 
NO2  concentration  schedules  are  very  similar  to  values  measured  outside  a research  house  in 
Chicago  (Figure  3.2  of  14).  Fine  particles  and  TVOCs  are  not  discussed  in  the  EPA  documents. 
The  ambient  fine  particle  concentration  was  chosen  based  on  the  average  of  reported  average 
measurements  for  four  US  cities  (Table  4 of  15).  The  TVOC  concentration  chosen  is  in  the 


10 


middle  of  the  reported  range  of  10  to  21 1 )ig/m^  measured  at  68  sites  in  the  US  (16).  The  fine 
particle  and  TVOC  concentrations  were  assumed  to  be  constant  throughout  the  day. 

In  addition  to  the  ambient  concentrations  in  Table  7 that  served  as  the  boundary  conditions  for 
the  indoor  sources,  elevated  levels  of  CO,  coarse  particles,  and  NO2  were  simulated  in  order  to 
evaluate  the  effect  of  the  lAQ  control  technologies  on  pollutants  brought  into  residences  jfrom  the 
outdoors.  These  elevated  pollutant  concentrations  were  selected  based  on  review  of  US  EPA  air 
quality  documents  (11,  12,  13)  and  were  specified  per  the  schedules  in  Table  8. 


Table  7 - Outdoor  pollutant  concentration  schedules 


Hour  of  day 

0-7 

7-9 

9-17 

17-19 

19-24 

CO  (ppm) 

1 

2 

1.5 

3 

1.5 

NO2  (ppm) 

0.02 

0.04 

0.02 

0.04 

0.02 

Fine  particles  (pg/m^) 

13 

13 

13 

13 

13 

TVOCs  (pg/m') 

100 

100 

100 

100 

100 

Table  8 - Elevated  outdoor  pollutant  concentrations  schedule 


Hour  of  day 

0-7 

7-9 

9-17 

17-19 

19-24 

CO  (ppm) 

4 

8 

7 

12 

6 

NO2  (ppm) 

0.2 

0.4 

0.2 

0.4 

0.2 

Coarse  particles  (|ig/m^) 

75 

75 

75 

75 

75 

The  Phase  I report  (2)  described  the  pollutant  sources  considered  for  inclusion  in  the  study.  The 
pollutant  sources  used  in  the  baseline  simulations  included  several  VOC  burst  sources  (medium 
strength  source  based  on  a polish  and  high  strength  source  based  on  a spray  carpet  cleanser  (17)), 
a constant  VOC  area  source  (based  on  a P VC  flooring  material  with  high  emissions  (18)),  and 
combustion  sources  (based  on  medium  source  strengths  for  ovens  and  space  heaters  (19))  of  CO, 
NO2,  and  fine  particles.  While  the  source  strength  used  for  the  flooring  material  isbased  on  a 
material  with  high  emissions,  it  is  only  moderately  higher  than  the  range  of  0.17  to  2.1 1 mg/m^h 
recently  reported  in  5 day  emission  tests  of  finished  particleboard  (20).  Table  9 lists  detailed 
information  on  these  sources  including  the  zones  (see  Figures  1 and  2 for  zone  labels,  also  BMT 
is  the  basement  zone)  in  which  they  are  located,  source  strengths,  and  schedules. 
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Table  9 - Pollutant  sources 


Source  name 

Pollutant 

Zone(s) 

Source  strength 

Schedule 

Burst  (medium) 

TVOCs 

Several 

300  mg/h 

9 - 9:30  am 

7 - 7:30  p.m. 

Burst  (high) 

TVOCs 

GAR  and  BMT 

1100  mg/h 

9 - 10  am 

7 - 8 p.m. 

Flooring  material 

TVOCs 

All  butGAR,ATC 

7.0  mg/h  m^ 

constant 

Oven 

CO 

KIT  (ranch  house), 

KFA  (two-story  house) 

1900  mg/h 

7 - 7:30  am 

6-7  p.m. 

Oven 

NO2 

KIT  (ranch  house), 

KFA  (two-story  house) 

160  mg/h 

7 - 7:30  am 

6-7  p.m. 

Oven 

Fine  particles 

KIT  (ranch  house), 

KFA  (two-story  house) 

0.2  mg/h 

7 -7:30  am 

6-7  p.m. 

Heater 

CO 

GAR  and  BMT 

1000  mg/h 

7 - 10  am  (GAR) 

7 - 9 p.m.  (BMT) 

Heater 

NO2 

GAR  and  BMT 

250  mg/h 

7 - 10  am  (GAR) 

7 - 9 p.m.  (BMT) 

Heater 

Fine  particles 

GAR  and  BMT 

2 mg/h 

7- 10  am  (GAR) 

7-9  p.m.  (BMT) 

In  addition  to  the  sources  listed  in  Table  9,  the  simulation  plan  in  the  Phase  I report  (2)  included 
a newly-finished  floor  as  a floor-area  based  decaying  source  of  VOCs.  A test  simulation  with  a 
medium  strength  source,  modeled  as  a first  order  exponential  decay  source  with  initial  emission 
rate  of  1 7400  mg/m^h  and  decay  constant  of  1 .24  h'^  (based  on  a stain  product  (21))  was 
performed  for  the  Miami  ranch  house.  This  source  resulted  in  extremely  high  concentrations  of 
TVOCs  with  a peak  concentration  of  over  2 g/m^  and  a concentration  of  37  mg/m^  at  the  end  of 
the  day.  None  of  the  lAQ  control  retrofits  being  evaluated  can  be  expected  to  have  a significant 
impact  on  the  extremely  high  concentrations  from  this  source  during  the  one-day  simulation 
period.  Therefore,  this  source  was  not  included  in  the  remaining  baseline  simulations.  Decaying 
high-strength  sources  such  as  this  one  are  of  interest  and  may  be  studied  in  the  future  with 
simulations  of  longer  duration. 

Reversible  sink  effects  for  the  VOCs  were  modeled  with  sink  elements  based  on  a boundaiy  layer 
diffusion  controlled  (BLDC)  model  with  a linear  adsorption  isotherm.  The  BLDC  adsorption 
model  is  described  by  Axley  (22).  The  parameters  required  for  this  sink  model  are  the  film  mass 
transfer  coefficient,  the  adsorbent  mass,  and  the  isotherm  partition  coefficient,  and  these 
parameters  would  vary  over  time  and  by  location  within  a house.  However,  since  little  real  data  is 
available  for  these  parameters  (which  depend  on  factors  such  as  gas  diffusion  properties,  airflow 
rates,  and  adsorbent  material)  and  because  the  goal  was  to  obtain  a reasonable  estimate  of  the 
reversible  sink  effects,  constant  values  were  used  for  all  of  the  parameters  and  only  the  adsorbent 
mass  was  varied  by  zone.  The  film  mass  transfer  coefficient  used  was  35  |xm/s  and  was 
calculated  from  equation  3.17a  of  Axley  (22)  with  an  assumed  air  velocity  of  0.001  m/s,  effective 
length  of  4 m,  Schmidt  number  of  1.0,  and  binary  diffusion  coefficient  of  1.0  x 10'^  mVs.  The 
partition  coefficient  used  was  0.5  g-air/g-sorbent  and  was  estimated  from  parameters  reported  for 
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an  empirical  sink  model  for  an  experimental  case  of  alkanes  emitted  by  a wood  stain  in  a test 
house  (23).  The  adsorbent  mass  used  was  based  on  a mass  of  6 kg  per  m^  of  adsorbent  surface 
area  which  was  assumed  to  be  equal  to  half  of  the  zone  interior  surface  area. 

Nitrogen  dioxide  decay  and  particle  deposition  were  modeled  as  single-reactant  first  order 
reactions  with  a single,  constant  value  in  all  rooms  of  the  houses.  Nitrogen  dioxide  decay 
depends  strongly  on  the  materials  present  in  a house  (e.g.,  floor  and  wall  coverings,  furnishings, 
etc.)  and  a wide  range  of  measured  values  have  been  reported  including  a range  of  0.09  - 13.74 
h*'  by  Lee  et  al.  (24).  Other  studies  have  reported  average  NO2  decay  rates  of  0.17,  0.29,  0.65, 

0.8,  0.82,  and  2.07  h*’  (25,  26,  27,  28, 29,  24).  The  kinetic  rate  coefficient  used  for  NO2  decay 
was  0.87  h'^  and  is  based  on  the  average  of  measurements  in  a contemporary  research  house 
reported  by  Leslie  et  al.  (14). 

Particle  deposition  depends  on  the  size  and  type  of  particles,  particle  concentration,  airflow 
conditions,  and  surfaces  available  for  deposition.  The  particle  decay  rate  used  for  fine  particles 
was  0.08  h‘^  and  was  reported  by  Traynor  et  al.  (30)  for  combustion  products  from  a 
wood-burning  stove  in  a test  house.  Offerman  et  al.  (31)  reported  a similar  mass-averaged  value 
of  0.1  h’^  for  tobacco  smoke  particles  in  a research  house.  The  decay  rate  can  be  calculated  as  the 
product  of  an  average  deposition  velocity  and  a room  surface-to-volume  ratio.  Assuming  a room 
surface-to-volume  ratio  of  2 m‘^  (the  actual  value  will  depend  on  room  geometry,  furnishings,  and 
surface  finishes),  a decay  rate  of  0.08  h‘^  corresponds  to  a deposition  velocity  of  approximately 
0.001  cm/s.  Sinclair  et  al.  (32,  33)  reported  higher  average  deposition  velocities  of  0.005  cm/s  for 
fine-mode  sulfate  in  telephone  equipment  buildings.  However,  the  nature  of  the  indoor 
environment,  and  especially  the  airflow  conditions,  in  a detached  single-family  home  and  a 
commercial  building  are  very  different.  Nazaroff  et  al.  (34)  discusses  the  use  of  deposition 
velocity  and  warns  that  "Deposition  velocities  determined  for  one  indoor  environment  can  only 
be  applied  to  another  to  the  extent  that  the  air  flow  conditions  are  similar." 

In  the  only  report  of  coarse  particle  deposition  rates  in  a test  house  found  in  the  literature,  Byrne 
et  al.  (35)  reported  values  of  1.51  and  2.10  h"'  for  4 pm  particles  in  an  unfurnished  and  furnished 
room,  respectively.  The  reported  mean  deposition  velocities  of  0.027  to  0.038  cm/s  fall  within 
the  range  of  approximately  .01  to  0.1  cm/s  calculated  from  a natural  convection  deposition  model 
by  Nazaroff  and  Cass  (36).  The  actual  decay  rate  for  the  coarse  outdoor  air  particles  modeled  in 
the  simulations  would  depend  on  the  size  distribution  of  the  particles.  Since  no  specific 
distribution  has  been  assumed , a decay  rate  of  1 .5  h'^  was  chosen  based  on  the  lower  value 
reported  by  Byrne. 

Baseline  Simulation  Results 

The  results  of  each  of  the  24  simulations  listed  in  Table  1 include  pollutant  concentrations  for  up 
to  18  pollutants  in  each  of  the  building  zones  for  each  15 -minute  time  step  of  the  24-hour 
simulation  period.  The  complete  transient  simulation  results  are  not  presented  here  but  are 
available  in  spreadsheet  files.  Instead,  this  section  presents  examples  of  the  transient  pollutant 
concentrations  for  selected  simulations.  Figures  3 through  6 show  the  pollutant  concentrations  in 
Zone  LDA  resulting  from  selected  pollutant  sources  for  simulation  SIMIFLC  (the  typical  Miami 
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ranch  house  in  cold  weather).  Figures  7 through  10  show  the  corresponding  results  for  SIMIFTC 
(the  tight  Miami  ranch  house  in  cold  weather).  Although  these  figures  are  only  examples  of 
transient  results,  some  observations  can  be  made  into  trends  and  factors  affecting  the  predicted 
contaminant  concentrations.  A complete  summary  of  peak  and  24-hour,  4-hour,  and  1 -horn- 
average  concentrations  for  all  baseline  simulations  are  included  in  Tables  la  through  24e  of 
Appendix  B. 

Figures  3 and  7 show  the  CO  concentrations  in  Zone  LDA  resulting  from  the  oven  and  heater 
sources  (CO.l  and  C0.2  in  Tables  1 and  4 of  Appendix  B)  for  SIMIFLC  and  SIMIFTC, 
respectively.  Both  graphs  show  two  daily  peaks  due  to  the  operation  of  the  oven.  For  the  tight 
house,  the  peaks  are  shifted  to  a slightly  later  time  due  to  reduced  outdoor  airflow  into  the  house 
which  resulted  in  less  mixing  of  CO  from  the  kitchen  into  the  rest  of  the  house.  The  heater 
source  causes  much  lower  concentrations  than  the  oven  source  due  to  the  low  airflow  rate  from 
the  garage  into  the  house  and  the  lower  source  emission  rate.  The  resulting  CO  concentrations  for 
the  heater  source  are  influenced  primarily  by  the  outdoor  level.  There  is  more  damping  of  the 
outdoor  variations  in  the  tight  house  than  the  typical  house  due  to  the  reduced  air  infiltration  rate. 

Figures  4 and  8 show  the  NO2  concentrations  (N02.1  and  N02.2  in  Tables  1 and  4 of  Appendix 
B)  in  Zone  LDA  resulting  from  the  oven  and  heater  sources  for  SIMIFLC  and  SIMIFTC.  The 
NO2  concentrations  show  some  of  the  same  characteristics  as  the  CO  concentrations,  with  two 
peaks  from  the  oven  and  a damped  dependence  on  the  outdoor  concentration  for  the  heater. 
However,  the  NO2  peaks  in  the  tight  house  are  significantly  less  than  in  the  typical  house  despite 
the  reduced  outdoor  airflow  into  the  house.  In  fact,  the  whole-house  average  NO2  concentrations 
are  lower  in  the  tight  house  than  in  the  typical  house  (0.025  ppm  vs.  0.026  ppm  for  the  oven 
source,  and  0.003  ppm  vs.  0.008  ppm  for  the  heater  source).  These  results  may  seem  surprising 
as  one  might  expect  the  reduced  infiltration  to  result  in  higher  NO2  concentrations  in  the  tight 
house.  However,  these  results  may  be  explained  by  the  impact  of  NO2  decay.  The  NO2  in  the 
house  is  either  generated  by  the  indoor  sources  or  brought  in  from  outside.  During  much  of  the 
day,  when  the  combustion  appliances  are  not  operated,  the  outdoor  air  is  the  source  of  indoor 
NO2  and,  due  to  NO2  decay,  the  indoor  concentrations  are  lower  than  the  outdoor  concentration. 
Therefore,  reducing  the  infiltration  actually  results  in  lower  indoor  NO2  concentrations.  When 
there  is  an  indoor  source  of  NO2,  a lower  infiltration  rate  may  still  result  in  lower  NO2 
concentrations  in  the  zones  without  the  source.  However,  the  source-zone  will  have  higher  NO2 
concentrations  when  the  infiltration  rate  is  lower  (the  peak  kitchen  concentration  from  the  oven  is 
1 .686  ppm  for  the  tight  house  and  1 .434  ppm  for  the  typical  house).  It  is  important  to  note  that 
this  result  of  lower  NO2  concentrations  in  tighter  houses  cannot  be  generalized  to  all  cases.  If  the 
NO2  decay  rate  was  lower,  the  indoor  NO2  generation  rate  was  higher,  or  the  outdoor  NO2 
concentration  was  lower,  the  tighter  house  could  have  higher  concentrations.  See  Tables  lb  and 
4b  of  Appendix  B for  the  peak  and  average  NO2  concentrations  results  for  these  simulations. 

The  TVOC  concentrations  in  Zone  LDA  resulting  from  the  burst  source  located  in  the  LDA  Zone 
(VOC4  in  Tables  1 and  4 of  Appendix  B)  are  shown  in  Figures  5 and  9 for  SIMIFLC  and 
SIMIFTC.  As  expected,  this  source  results  in  two  daily  peaks  due  to  the  source  schedule  and 
higher  concentrations  in  the  tight  house  due  to  the  reduced  airflow  into  the  house. 
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Figures  6 and  10  show  the  CO  and  coarse  particle  concentrations  (C0.3  and  PART.3  of  Tables  1 
and  4 in  Appendix  B)  for  SEMIFLC  and  SEMI  FTC,  respectively,  due  to  the  elevated  outdoor 
pollutant  concentrations  of  Table  8.  The  CO  concentration  in  the  typical  house  tracks  the  outdoor 
concentration  with  a time  lag  based  on  the  building  time  constant  (related  to  the  inverse  of  the 
building  air  change  rate).  The  particle  concentration  shows  the  effect  of  the  ETVAC  system 
cycling  which  changes  the  air  change  rate  of  the  house  and  filters  particles  from  the  air.  When  the 
furnace  is  on,  the  concentration  of  coarse  particles  decreases  due  to  the  impact  of  the  furnace 
filter.  The  tight  house  results  exhibit  damped  CO  peaks  and  valleys  due  to  the  longer  time 
constants.  Because  the  particles  come  from  outdoors,  the  lower  air  change  rates  result  in  lower 
particle  concentrations. 
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Figure  3 - CO  in  Zone  LDA  for  SiM1  FLC 
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Rgure  5 - TVOCs  in  Zone  LDA  for  SIM1 FLC 
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Figure  7 - CO  in  Zone  LDA  for  SIM1  FTC 
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Rgure  8 - N02  in  Zone  LDA  for  SIM1  FTC 
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Rgure  9 - TVOCs  in  Zone  LDA  for  SIM1  FTC 
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Indoor  Air  Quality  Controls 

This  section  describes  the  indoor  air  quality  control  technologies  that  will  be  evaluated  in  the 
study.  These  technologies  will  be  incorporated  into  the  baseline  house  models  to  determine  their 
effectiveness  in  controlling  the  selected  pollutant  sources.  The  three  technologies  described  in 
this  section  include  the  following: 

Electrostatic  particulate  filtration 
Heat  recovery  ventilation 

Outdoor  air  intake  damper  on  the  forced-air  system  return 

This  section  describes  each  of  these  technologies  and  includes  revisions  of  the  baseline  house 
duct  drawings.  In  addition,  this  section  contains  an  estimate  of  the  equipment  and  installation 
costs  and  a revision  of  the  thermal  load  calculations  based  on  the  modifications.  Finally,  the 
impacts  of  each  of  these  technologies  on  “other  contaminants”  are  discussed.  These  other 
contaminants,  as  described  in  the  original  project  work  statement,  include  contaminants  that  have 
typically  been  of  concern  to  designers  of  residential  ventilation  systems  including  cooking  odors, 
tobacco  smoke,  moisture,  outdoor  pollen,  outdoor  odors,  and  ozone. 

Electrostatic  Particulate  Filtration 

The  first  lAQ  control  technology  is  increased  particulate  filtration  through  the  installation  of 
passive,  electrostatic  particulate  filters.  These  filters  were  chosen  based  on  the  availability  of 
performance  data.  In  addition,  the  low  pressure  drop  through  these  filters  enables  their 
installation  without  modification  of  the  existing  forced-air  distribution  system.  The  baseline 
houses  are  assumed  to  have  standard  furnace  filters  with  an  ASHRAE  dust  spot  efficiency  of  less 
than  20%  and  an  arrestance  of  90%.  These  values  are  based  on  tests  conducted  in  accordance 
with  ASHRAE  Standard  52.1  (37).  The  increased  filtration  is  based  on  the  use  of  electrostatic 
filters  with  an  ASHRAE  dust  spot  efficiency  of  30%  and  an  arrestance  of  95%. 

Although  the  efficiencies  of  particulate  filters  change  over  time  as  they  become  loaded,  the 
computer  simulations  in  this  project  will  employ  a constant  filter  efficiency.  The  efficiencies  of 
the  baseline  and  improved  filters  used  in  the  simulations  will  be  as  follows: 


Baseline 

Control  #1 

Particles  <2.5  pm  in  diameter 

5% 

30% 

Particles  between  2.5  and  10  pm  in  diameter 

90% 

95% 

The  improved  filters  are  installed  in  place  of  the  regular  furnace  filters.  Their  location  is 
indicated  in  the  revised  duct  drawings  showing  all  of  the  LAQ  control  technologies.  Figures  13 
through  16. 
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The  installation  of  the  improved  filters  are  assumed  not  to  affect  the  thermal  loads  of  the  houses. 
Due  to  a higher  pressure  drop  through  the  filters,  they  may  cause  a slight  reduction  in  the  airflow 
rate  through  the  system,  which  could  affect  the  pressures  across  the  building  envelope  and  the 
resultant  building  infiltration  rates.  However,  this  effect  is  expected  to  be  small,  and  the  thermal 
load  calculations  were  not  modified  for  this  control  technology. 

The  cost  of  this  first  control  technology  includes  the  cost  of  the  filters  themselves  and  their 
installation.  For  comparison,  the  furnace  filters  in  the  baseline  houses  are  assumed  to  cost  $2 
each  and  to  be  changed  every  month.  Therefore,  the  annual  cost  of  the  baseline  filters  is  $24.  The 
improved  filters  are  assumed  to  cost  $15  each  and  to  be  changed  every  2 months.  Therefore,  the 
annual  cost  of  the  improved  filters  is  $90. 

The  installation  of  improved  filters  will  reduce  the  concentrations  of  the  so-called  “other 
contaminants”  in  the  houses  to  the  degree  that  the  filtration  of  each  contaminant  is  increased.  The 
concentrations  of  particulate  contaminants  with  outdoor  sources  (pollen)  will  be  reduced  due  to 
the  increased  particulate  filtration.  The  concentrations  of  VOCs  associated  with  outdoor  odors 
will  not  be  decreased.  The  increased  filtration  will  not  affect  indoor  ozone  levels  due  to  outdoor 
sources,  since  ozone  removal  rates  will  be  unaffected  by  the  new  filters.  In  addition,  these 
electrostatic  filters  are  not  sources  of  ozone  themselves.  The  concentrations  of  other 
contaminants  with  indoor  sources  will  also  be  affected  to  the  degree  that  the  filtration  of  each 
contaminant  is  increased.  The  levels  of  cooking  odors  and  tobacco  smoke  will  be  decreased 
based  on  the  increased  filter  efficiency  for  both  fine  and  coarse  particulates.  Indoor  moisture 
levels  will  be  unaffected  by  the  new  filters  because  the  outdoor  air  change  rates  will  not  be 
affected  and  because  the  improved  filters  have  no  humidification  or  dehumidification  impacts. 

Electronic  air  cleaners  are  also  of  interest  and  may  be  investigated  in  follow-up  work.  The 
existence  of  reliable  performance  data  is  being  investigated. 

Heat  Recovery  Ventilator 

The  second  lAQ  control  technology  is  the  installation  of  a heat  recovery  ventilator  (HRV)  in 
conjunction  with  the  forced-air  distribution  system.  As  seen  in  Figure  1 1,  the  device  brings 
outdoor  air  into  the  building  where  it  exchanges  heat  with  an  airstream  leaving  the  return  side  of 
the  forced  air  system.  Under  heating  conditions,  the  outdoor  air  is  warmed  by  the  outgoing 
airstream,  and  under  cooling  the  outdoor  air  is  cooled.  The  outgoing  airstream  is  exhausted  to  the 
outdoors  after  leaving  the  heat  recovery  ventilator.  The  airstream  from  outdoors  flows  into  the 
return  side  of  the  forced-air  system  after  leaving  the  HRV. 
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Figure  1 1 - Schematic  of  Heat  Recovery  Ventilator 

The  HRV  specifications  are  based  on  a commercially-available  model  designed  for  residential 
use  and  installation  in  conjunction  with  forced-air  systems.  The  airflow  rate  capacity  of  the 
device  was  selected  to  obtain  an  air  change  rate  of  at  least  0.5  air  changes  per  hour  (ach)  at  full 
flow.  The  actual  outdoor  airflow  rate  during  operation  was  selected  to  provide  0.35  ach  through 
the  HRV.  The  actual  whole  building  air  change  rate  will  also  include  envelope  infiltration,  which 
in  turn  depends  on  the  airtightness  of  the  house,  weather  conditions  and  ventilation  equipment 
operation.  The  HRV  specifications  for  the  four  houses  are  as  follows: 

Miami,  2-story 

Airflow  capacity:  30  to  60  L/s  (65  to  127  cfm),  roughly  0.25  to  0.5  ach 
Airflow  rate  during  operation:  44  L/s  (93  cfin) 

Efficiency:  69%  at  0 °C  (32  °F),  60%  at  -25  °C  (-13  °F) 

Maximum  power  consumption:  1 1 5 W 
No  defrost 

Miami,  Ranch 

Airflow  capacity:  30  to  60  L/s  (65  to  127  cfrn),  roughly  0.4  to  0.8  ach 
Airflow  rate  during  operation:  26  L/s  (55  cfrn) 

Efficiency:  69%  at  0 °C  (32  °F),  60%  at  -25  °C  (-13  °F) 

Maximum  power  consumption:  1 1 5 W 
No  defrost 

Minneapolis,  2-story 

Airflow  capacity:  55  to  95  L/s  (1 15  to  200  cfrn),  roughly  0.3  to  0.5  ach 
Airflow  rate  during  operation:  66  L/s  (140  cfin) 

Efficiency:  68%  at  0 °C  (32  °F),  61%  at  -25  °C  (-13  °F) 

Maximum  power  consumption:  216  W 
Defrost  cycle 
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Minneapolis,  Ranch 

Airflow  capacity:  30  to  70  L/s  (65  to  150  cfin),  roughly  0.2  to  0.5  ach 
Airflow  rate  during  operation:  52  L/s  (1 10  cfin) 

Efficiency:  76%  at  0 °C  (32  °F),  56%  at  -25  °C  (-13  °F) 

Maximum  power  consumption:  105  W 
Defi-ost  cycle 

The  defi-ost  cycle  involves  closing  the  outdoor  air  dampers  for  5 minutes  when  the  outdoor 
temperature  is  below  -5  °C  (23  °F).  For  outdoor  temperatures  between  -5  and  -30  °C  (23  and  -22 
°F),  each  5-minute  defrost  cycle  is  followed  by  a 35  minute  period  of  air  exchange  before  the 
next  defi-ost  cycle.  For  outdoor  temperatures  below  -30  °C  (-22  °F),  each  5-minute  defrost  cycle 
is  followed  by  20  minutes  of  air  exchange. 

The  HRV  can  be  operated  in  several  different  control  modes.  The  operation  of  the  device  and  the 
fan  speed  (high  or  low)  can  be  controlled  by  a timer,  manually  by  the  occupant  or  by  a 
dehumidistat. 

The  installation  of  the  HRV  in  each  of  the  four  houses  is  indicated  in  the  revised  duct  drawings 
in  Figures  13  through  16. 

The  thermal  loads  of  the  houses  are  affected  by  the  installation  and  operation  of  the  HRV  due  to 
the  increased  outdoor  air  change  rate  of  the  house  when  the  devices  are  in  operation.  The  air 
change  rate  due  to  the  HRV  operation  is  assumed  to  be  additive  to  the  baseline  infiltration  rate  of 
0.75  ach  assumed  for  the  design  thermal  load  calculations.  The  thermal  loads  are  increased  by 
only  a fi-action  of  the  increased  outdoor  air  change  rate  based  on  the  heat  exchange  efficiencies  of 
the  devices.  For  an  additional  air  change  rate  of  0.35  ach  and  the  rated  heat  exchange  efficiencies 
of  the  HRVs,  the  revised  design  thermal  loads  for  the  four  houses  are  given  below.  The  baseline 
design  thermal  loads  are  described  in  detail  in  the  Phase  I report  (2). 


Miami,  2-story 

Baseline 

With  HRV 

Heating 

2.87  kW 

3.14  kW 

Cooling 

6.43  kW 

6.60  kW 

Miami,  Ranch 

Baseline 

With  HRV 

Heating 

1.83  kW 

1.99  kW 

Cooling 

5.76  kW 

5.88  kW 

Minneapolis,  2-story 

Baseline 

With  HRV 

Heating 

12.64  kW 

13.59  kW 

Cooling 

6.21  kW 

6.36  kW 

Minneapolis,  Ranch 

Baseline 

With  HRV 

Heating 

9.25  kW 

9.86  kW 

Cooling 

4.89  kW 

4.97  kW 
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The  cost  of  the  HRVs  includes  the  cost  of  the  equipment  and  installation,  the  operating  costs  for 
the  fans  in  the  devices  and  the  increased  energy  consumption  due  to  the  additional  outdoor  air 
change  of  the  building.  The  cost  of  the  equipment  is  $500  for  both  of  the  Miami  houses,  $600  for 
the  Minneapolis  ranch  house  and  $700  for  the  Minneapolis  two-story  house.  These  are  list  prices 
from  the  manufacturer  of  the  HRV  on  which  the  specifications  are  based.  The  installation  costs 
are  more  variable,  based  on  the  layout  of  the  house  and  local  labor  rates,  and  they  can  range  from 
$200  to  $500.  The  cost  of  the  energy  consumed  by  the  device  and  by  the  additional  outdoor  air 
change  rate  requires  detailed  thermal  modeling  of  the  building  and  system.  As  discussed  in  the 
Phase  I report  of  the  project,  such  modeling  is  beyond  the  scope  of  this  project. 

The  installation  of  the  HRV  will  impact  the  so-called  “other  contaminants”  in  the  houses  due  to 
the  increased  outdoor  air  change  rate.  Due  to  the  additional  outdoor  airflow  into  the  houses,  the 
concentrations  of  contaminants  with  outdoor  sources  (pollen,  outdoor  odors  and  ozone)  will 
increase.  For  a simple,  nonreactive  and  unfiltered  contaminant,  there  will  be  an  increased 
contaminant  load  equal  to  the  outdoor  concentration  multiplied  by  the  outdoor  airflow  rate.  The 
impact  of  particulates  will  be  reduced  based  on  the  efficiency  of  the  filters  in  the  HRV  and  of  the 
furnace  filter.  The  impact  of  outdoor  ozone  will  be  reduced  somewhat  by  losses  on  the  interior 
surfaces  of  the  HRV  ductwork.  The  concentrations  of  other  contaminants  with  indoor  sources 
(cooking  odors  and  tobacco  smoke)  will  be  reduced  based  on  the  increased  air  change  rate  of  the 
building.  The  impact  of  the  additional  ventilation  on  moisture  will  depend  on  the  building 
location,  indoor  moisture  sources,  and  season.  Indoor  humidity  levels  will  be  reduced  when  there 
are  large  indoor  sources  and  low  relative  humidity  outdoors,  but  will  be  increased  when  the 
outdoor  humidity  is  higher  than  the  indoor  level.  Detailed  modeling  of  moisture  transport  is 
required  to  assess  these  impacts  and  is  beyond  the  scope  of  the  current  project. 

Outdoor  Intake  Duct 

The  third  lAQ  control  technology  is  the  installation  of  an  outdoor  air  intake  duct  on  the  return 
side  of  the  forced  air  distribution  system.  As  seen  in  Figure  12,  the  system  consists  of  an  intake,  a 
duct,  a motorized  damper,  and  a volume  damper  for  adjusting  the  airflow  rate,  and  is  connected 
to  the  return  side  of  the  return  duct.  The  maximum  airflow  rate  capacity  of  the  intake  is  78  L/s 
(165  cfin),  which  corresponds  to  the  following  air  change  rates  for  the  four  houses: 

Miami,  2-story:  0.62  ach 

Miami,  Ranch:  1 .05  ach 

Minneapolis,  2-story:  0.41  ach 

Minneapolis,  Ranch:  0.53  ach 
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Figure  12  - Schematic  of  Outdoor  Air  Intake  Duct 

The  actual  airflow  rate  through  the  intake  depends  on  the  position  of  the  volume  damper,  the 
overall  airflow  resistance  of  the  intake  system,  and  the  pressure  developed  by  the  forced-air  fan. 
In  the  computer  simulations,  it  is  assumed  that  the  volume  damper  is  adjusted  such  that  the 
intake  system  provides  0.35  ach  to  the  building  when  the  furnace  fan  is  in  operation.  This  air 
change  rate  corresponds  to  the  following  outdoor  air  intake  rates  for  the  four  buildings: 

Miami,  2-story:  44  L/s  (93  cfin) 

Miami,  Ranch:  26  L/s  (55  cfm) 

Minneapolis,  2-story:  66  L/s  (140  cfm) 

Minneapolis,  Ranch:  52  L/s  (1 10  cfin) 

The  motorized  damper  can  be  controlled  in  several  different  ways.  It  is  generally  interlocked  with 
the  forced-air  system  fan  so  that  it  opens  only  when  the  forced-air  fan  is  operating.  The 
motorized  damper  can  also  be  controlled  to  open  based  on  a timer,  dehumidistat  or  pollutant  (e.g. 
carbon  monoxide  or  carbon  dioxide)  sensor. 

The  installation  of  the  outdoor  air  intake  duct  in  each  of  the  four  houses  is  indicated  in  the 
revised  duct  drawings  in  Figures  13  through  16. 

The  thermal  loads  of  the  houses  are  affected  by  the  installation  and  operation  of  the  outdoor  air 
intake  duct  due  to  the  increased  outdoor  air  change  rate  of  the  house  when  the  devices  are  in 
operation.  The  air  change  rate  due  to  the  HRV  operation  is  assumed  to  be  additive  to  the  baseline 
infiltration  rate  of  0.75  ach  assumed  for  the  design  thermal  load  calculations.  Based  on  an 
additional  air  change  rate  of  0.35  ach  and  no  heat  exchange,  the  design  thermal  loads  for  the  four 
houses  are  given  below.  The  baseline  thermal  loads  were  described  in  detail  in  the  Phase  I report 
(2). 
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Miami,  2-story 

Baseline 

With  OAID 

Heating 

2.87  kW 

3.54  kW 

Cooling 

6.43  kW 

6.96  kW 

Miami,  Ranch 

Baseline 

With  OAID 

Heating 

1.83  kW 

2.23  kW 

Cooling 

5.76  kW 

6.09  kW 

Minneapolis,  2-story 

Baseline 

With  OAID 

Heating 

12.64  kW 

15.00  kW 

Cooling 

6.21  kW 

6.71  kW 

Minneapolis,  Ranch 

Baseline 

With  OAID 

Heating 

9.25  kW 

10.73  kW 

Cooling 

4.89  kW 

5.18  kW 

The  cost  of  the  outdoor  air  intake  duct  includes  the  cost  of  the  equipment  and  installation  and  the 
increased  energy  consumption  due  to  the  additional  outdoor  air  change  of  the  building.  The  cost 
of  the  equipment,  including  the  contiols  and  the  motorized  dampers,  is  $750  based  on  list  prices 
from  the  manufacturer  of  the  outdoor  air  intake  duct  on  which  the  specifications  are  based.  The 
installation  costs  are  more  variable,  based  on  the  layout  of  the  house  and  local  labor  rates,  and 
they  can  range  from  $100  to  $300.  The  cost  of  the  energy  consumed  by  the  device  and  by  the 
additional  outdoor  air  change  rate  requires  detailed  thermal  modeling  of  the  building  and  system. 
As  discussed  in  the  Phase  I report  of  the  project  (2),  such  modeling  is  beyond  the  scope  of  this 
project. 

The  installation  of  the  outdoor  air  intake  duct  will  impact  the  so-called  “other  contaminants”  in 
the  houses.  Due  to  the  additional  outdoor  airflow  into  the  houses,  the  concentrations  of 
contaminants  with  outdoor  sources  (pollen,  outdoor  odors  and  ozone)  will  increase.  For  a simple, 
nonreactive  and  unfiltered  contaminant,  the  impact  will  be  an  increased  contaminant  load  equal 
to  the  outdoor  concentration  multiplied  by  the  outdoor  airflow  rate.  The  impact  of  particulates 
will  be  lessened  based  on  the  removal  efficiency  of  the  furnace  filter.  The  impact  of  ozone  will 
be  lessened  by  losses  on  the  interior  surfaces  of  the  ductwork.  The  concentrations  of  other 
contaminants  with  indoor  sources  (cooking  odors  and  tobacco  smoke)  will  be  reduced  based  on 
the  increased  air  change  rate  of  the  building.  As  in  the  case  of  the  HRV,  the  impact  of  the 
additional  ventilation  on  moisture  will  depend  on  the  building  location,  indoor  moisture  sources, 
and  season. 
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Figure  13  - lAQ  Controls  for  Minneapolis  Ranch  House 
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Figure  14  - lAQ  Controls  for  Minneapolis  2-Story  House 


HRV  HRV 


Inti 

ce  Exhi 

ust 

1 

SUEEU^TRUNIC 

1 

r ^ CONTROL#! 

\ HEAT  RECOVERY 

\ VENTILArOR 

CONTROL  #3 

OUTDOOR  AIR 

INTAKE 

ATTIC 


Figure  15  - lAQ  Controls  for  Miami  Ranch  House 
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Figure  16  - lAQ  controls  for  Miami  2-Story  House 


Preliminary  Simulation  of  lAQ  Control  Retrofits 

This  section  describes  the  preliminary  simulations  of  the  lAQ  control  retrofits.  These  simulations 
involved  modifying  selected  baseline  simulation  cases  with  the  lAQ  control  retrofits  described 
above.  The  preliminary  simulations  were  performed  to  verify  the  ability  of  the  program  to  model 
the  control  technologies.  In  Phase  II.B  of  the  study,  all  of  the  baseline  cases  tvill  be  modified  to 
incorporate  each  of  the  lAQ  control  retrofits. 

lAQ  Control  Retrofits 

The  lAQ  control  retrofits  selected  for  modeling  in  this  study  are  an  electrostatic  particulate  filter, 
a heat  recovery  ventilator,  and  an  outdoor  air  intake  damper  installed  on  the  forced-air  system 
return.  These  three  technologies  were  described  in  detail  in  the  previous  section.  This  section 
discusses  only  the  details  important  to  modeling  them  with  CONTAM93  (1). 

The  electrostatic  particulate  filters  selected  for  the  study  have  a filter  efficiency  of  30%  for  fine 
particles  (emitted  by  the  combustion  sources  in  these  simulations)  and  95%  for  coarse  particles 
(associated  with  the  elevated  outdoor  air  concentrations).  The  filters  will  be  modeled  by 
replacing  the  standard  furnace  filters  in  the  baseline  HVAC  systems  with  the  electrostatic  filters. 
The  filter  efficiency  will  be  modeled  as  constant  over  time  and  impacts  on  airflow  through  the 
system  will  be  neglected. 

The  second  lAQ  control  retrofit  is  the  installation  of  a heat  recovery  ventilator  (HRV)  in 
conjunction  with  the  HVAC  system.  The  HRV  draws  air  from  the  return  side  of  the  forced-air 
system  and  replaces  it  with  outdoor  air  drawn  through  the  heat  exchanger.  The  outdoor  airflow 
rate  supplied  will  be  44  L/s  for  the  Miami  2-story  house,  26  L/s  for  the  Miami  ranch  house,  66 
L/s  for  the  Minneapolis  2-story  house,  and  52  L/s  for  the  Minneapolis  ranch  house.  The  HRV 
will  be  modeled  by  setting  the  outdoor  airflow  rate  for  each  HVAC  system  to  the  appropriate 
fraction  of  the  total  system  supply  airflow  rate.  Thus,  the  desired  amount  of  outdoor  air  will  be 
supplied  whenever  the  HVAC  system  is  operating.  The  HVAC  systems  will  be  operated  on  the 
same  schedules  determined  for  the  baseline  simulations  based  on  thermal  loads.  Other  possible 
control  options  (such  as  constant  operation  or  demand  control)  will  not  be  studied. 

Other  considerations  in  modeling  the  HRV  include  filtration  of  the  incoming  outdoor  air  and  the 
HRV  defrost  cycle.  A standard  furnace  filter  (with  efficiencies  of  5%  for  fine  particles  and  90% 
for  coarse  particles)  will  be  included  in  the  outdoor  air  intake  path  of  the  HRV.  The  HRV 
employs  a defrost  cycle  in  cold  weather  which  involves  periodically  closing  the  outdoor  air 
damper.  However,  operation  of  the  defrost  cycle  will  be  neglected  in  the  simulations. 

The  third  lAQ  control  retrofit  is  the  installation  of  an  outdoor  air  intake  duct  on  the  return  side  of 
the  HVAC  system,  which  draws  outdoor  air  into  the  return  side  of  the  forced-air  system 
whenever  it  is  operating.  This  retrofit  will  be  modeled  similar  to  the  HRV.  The  baseline  HVAC 
system  will  be  modified  to  include  a constant  fraction  of  outdoor  air  whenever  the  HVAC  system 
is  operating.  The  outdoor  air  supply  airflow  rates  will  be  the  same  as  listed  above  for  the  HRV, 
and  a standard  furnace  filter  will  be  included  in  the  outdoor  air  intake  path.  The  primary 
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difference  between  the  outdoor  air  intake  damper  and  the  HRV  is  that  the  outdoor  air  intake 
damper  does  not  include  an  exhaust  duct.  Therefore,  the  outdoor  airflow  will  tend  to  pressurize 
the  house.  This  effect  will  be  modeled  by  reducing  the  HVAC  return  flows  by  an  amount  equal  to 
the  outdoor  air  supplied  to  the  system. 

Results  of  Preliminary  Simulation  of  lAQ  Control  Retrofits 

The  baseline  case  selected  for  modification  with  the  LAQ  control  retrofits  was  SIMIFLC,  the 
Miami  ranch  house  with  typical  airtightness  in  cold  weather.  The  simulations  with  the 
electrostatic  particulate  filtration,  the  HRV,  and  the  outdoor  air  intake  damper  are  referred  to  as 
SIMIFLCF,  SIMIFLCH,  and  SIMIFLCO,  respectively. 

The  results  of  each  simulation  includes  pollutant  concentrations  for  up  to  18  pollutants  in  each  of 
the  building  zones  for  each  15  minute  time  step  of  the  24  hour  simulation  period.  As  was  the  case 
for  the  baseline  simulations,  the  complete  transient  simulation  results  are  not  presented  here  but 
are  available  in  spreadsheet  files.  Figures  17  through  20  show  examples  of  the  transient  pollutant 
concentrations  and  Tables  25a  through  27e  of  Appendix  B present  a complete  summary  of  peak 
and  24-hour,  4-hour,  and  1-hour  average  concentrations  for  the  preliminary  lAQ  control  retrofit 
simulations  as  described  for  the  baseline  simulations. 

Figure  17  shows  total  volatile  organic  compound  (TVOC)  concentrations  in  Zone  LDA  resulting 
from  the  constant  floor  source  (VOC2  of  Tables  1, 26,  and  27  of  Appendix  B)  for  SIMIFLC, 
SIMIFLCH,  and  SIMIFLCO.  Since  SIMIFLCF  differs  from  SIMIFLC  by  improved  particle 
filtration  efficiency,  all  VOC  concentrations  in  SIMIFLCF  are  identical  to  SIMIFLC  and  are  not 
shown.  Both  outdoor  air  intake  devices  result  in  modest  reductions  in  the  TVOC  concentrations 
in  the  zone,  with  the  HRV  having  a slightly  greater  effect.  The  HRV  may  have  a greater  effect 
because  it  has  a neutral  effect  on  indoor  pressure  (compared  to  the  outdoor  air  intake  damper 
which  pressurizes  the  building)  resulting  in  a greater  average  air  change  rate.  The  24-hour 
average  TVOC  concentration  in  Zone  LDA  is  6040  jxg/m^  5545  |ig/m^  and  5720  jig/m^  for 
SIMIFLC,  SIMIFLCH,  and  SIMIFLCO,  respectively  (see  Tables  Ic,  26c,  and  27c  of  Appendix 
B). 

Figure  18  shows  TVOC  concentrations  in  Zone  LDA  resulting  from  a burst  VOC  source  in  the 
garage  (VOC5  of  Tables  1, 26,  and  27  of  Appendix  B).  The  VOC  concentrations  for  SIMIFLC 
and  SIMIFLCH  are  nearly  identical  while  the  concentrations  for  SIMIFLCO  are  somewhat 
lower.  The  24  hour  average  TVOC  concentration  in  Zone  LDA  for  this  source  is  141  |Lig/m^,  140 
|Lig/m^,  and  132  for  SIMIFLC,  SIMIFLCH,  and  SIMIFLCO,  respectively  (see  Tables  Ic, 
26c,  and  27c  of  Appendix  B).  The  slightly  reduced  concentrations  for  SIMIFLCO  is  due  to  the 
effect  of  the  outdoor  air  pressurizing  the  interior  of  the  house  which  reduces  the  transport  of  the 
contaminant  from  the  garage. 

Figure  19  shows  fine  particle  concentrations  in  the  kitchen  resulting  from  the  oven  source 
(PART.l  of  Tables  1,  25,  26,  and  27  of  Appendix  B).  The  improved  filtration  in  case  SIMIFLCF 
resulted  in  lower  concentrations  while  the  outdoor  air  intake  devices  had  very  little  impact, 
possibly  because  the  outdoor  air  particle  concentration  of  13  |ig/m^  is  close  to  the  24  hour 
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average  baseline  concentration.  The  24  hour  average  fine  particle  concentration  in  Zone  KIT  is 
1 1.4  lLlg/m^  9.8  ^lg/m^  1 1.6  ^g/m^  and  1 1.6  jig/m'  for  SIMIFLC,  SIMIFLCF,  SIMIFLCH,  and 
SIMIFLCO,  respectively  (see  Tables  Ic,  25c,  26c,  and  27c  of  Appendix  B). 

Figure  20  shows  coarse  particle  concentrations  in  Zone  LDA  resulting  from  elevated  outdoor 
levels  (PART.3  of  Tables  1, 25,  26,  and  27  of  Appendix  B).  None  of  the  lAQ  control  retrofits 
resulted  in  a significant  impact  on  the  coarse  particle  concentrations.  The  24  hour  average  coarse 
particle  concentration  in  Zone  LDA  is  13.7  jig/m^  13.6  |Xg/m^  13.8  lLlg/m^  and  13.5  jig/m^  for 
SIMIFLC,  SIMIFLCF,  SIMIFLCH,  and  SIMIFLCO,  respectively  (see  Tables  Ic,  25c,  26c,  and 
27c  of  Appendix  B).  Possible  explanations  for  the  small  changes  include  the  relatively  small 
increase  in  filtration  efficiency  for  the  electrostatic  particulate  filter  (from  90%  to  95%)  and  the 
inclusion  of  a standard  filter  in  the  outdoor  air  intake  path  for  both  the  HRV  and  the  outdoor  air 
intake  damper.  The  outdoor  air  intake  filter  limits  the  number  of  particles  brought  in  with  the 
outdoor  air. 
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Figure  1 7 - TVOC  concentrations  in  zone  LDA  from  floor  source 
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Rgure  19  - Rne  particle  concentrations  in  zone  KIT  from  oven]| 
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Rgure  20  - Coarse  particle  concentrations  in  LDA  from  elevated  outdoor  levels 
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Summary 


The  National  Institute  of  Standards  and  Technology  (NIST)  has  completed  Phase  11.  A of  a project 
for  the  U.S.  Consumer  Product  Safety  Commission  (CPSC)  to  study  the  impact  of  HVAC 
systems  on  residential  indoor  air  quality  and  to  assess  the  potential  for  using  residential  forced-air 
systems  to  control  indoor  pollutant  levels.  In  this  effort,  NIST  is  performing  whole  building 
airflow  and  contaminant  dispersal  computer  simulations  with  the  program  CONTAM93  to  assess 
the  ability  of  modifications  of  central  forced-air  heating  and  cooling  systems  to  control  pollutant 
sources  relevant  to  the  residential  environment.  During  Phase  n.A  of  this  project,  three  major 
efforts  were  completed:  baseline  simulations  of  contaminant  levels  without  lAQ  controls,  design 
of  the  LAQ  control  retrofits,  and  preliminary  simulations  of  contaminant  levels  with  the  lAQ 
control  retrofits. 

It  is  important  to  note  that  the  project  is  essentially  a scoping  study  to  conduct  a preliminary 
assessment,  using  computer  simulation,  of  the  potential  for  using  forced-air  HVAC  systems  to 
improve  residential  lAQ.  The  project  results  are  also  limited  by  the  lack  of  high  quality  input  data 
for  some  simulation  inputs  and  the  lack  of  a thorough  empirical  evaluation  of  the  model's 
predictive  capability.  Despite  these  limitations,  the  project  is  expected  to  identify  key  issues  for 
further  analysis  and  experimental  work  to  meet  the  overall  goal  of  cost-effective  lAQ  control  in 
residential  buildings. 

This  report  described  the  input  data  used  to  model  the  baseline  houses  with  CONTAM93 
including  the  configuration  of  the  building  zones,  the  air  leakage  of  the  building  envelopes  and  of 
interior  partitions,  wind  pressure  profile  on  the  building  envelope,  pollutant  source  strengths  and 
temporal  profiles,  heating  and  cooling  system  flows,  furnace  filter  efficiency,  pollutant  sinks, 
pollutant  decay  or  deposition,  and  ambient  weather  and  pollutant  concentrations.  The  results  of 
the  baseline  simulations  including  transient  pollutant  concentrations  for  selected  simulations  and 
a summary  of  peak  and  average  concentrations  for  all  baseline  simulations  were  also  presented.  It 
should  be  noted  that  the  results  for  any  one  simulation  may  be  counter-intuitive  and  should  not 
be  generalized  to  all  cases. 

Three  indoor  air  quality  control  technologies  were  then  selected  for  incorporation  into  the 
baseline  house  models  to  determine  their  effectiveness  in  controlling  the  modeled  pollutant 
sources.  The  technologies  selected  include  the  following:  an  electrostatic  particulate  filter  with 
efficiencies  of  30%  for  fine  particles  and  95%  for  coarse  particles,  a heat  recovery  ventilator 
(HRV)  providing  an  actual  outdoor  airflow  of  0.35  ach,  and  an  outdoor  air  intake  damper  on  the 
forced-air  system  return  also  providing  an  actual  outdoor  airflow  of  0.35  ach.  The  annual  cost  of 
the  filters  was  estimated  at  $90.  The  estimated  installation  and  equipment  costs  of  the  HRV  and 
of  the  outdoor  air  intake  duct  were  $700  to  $1200  and  $850  to  $1050,  respectively.  Detailed 
thermal  modeling  of  the  building  and  system  would  be  required  to  determine  the  annual  energy 
costs  of  these  devices  and  is  beyond  the  scope  of  this  project. 

Selected  baseline  cases  were  then  modified  to  implement  these  LAQ  control  retrofits  and 
preliminary  simulations  were  performed  to  verify  the  ability  of  the  program  to  model  the  control 
technologies.  The  results  for  the  lAQ  control  retrofits  are  presented  as  examples  only  and  are  not 
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intended  to  be  used  to  evaluate  the  effectiveness  of  the  controls.  In  Phase  n.B  of  the  study,  all  of 
the  baseline  cases  will  be  modified  to  incorporate  each  of  the  lAQ  control  retrofits.  The  Phase 
n.B  simulation  results  will  be  compared  to  the  baseline  simulation  results  to  determine  the 
effectiveness  of  the  lAQ  control  technologies  at  reducing  contaminant  levels  in  single-family 
residential  buildings. 
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Appendix  A Airflow  Modeling  Results 

CONTAM93  was  used  to  analyze  airflow  in  the  houses  using  two  approaches:  simulated  fan 
pressurization  tests  and  directly  calculated  whole  building  air  change  rates  under  a range  of  wind 
speed  and  indoor  - outdoor  temperature  differences. 

Fan  pressurization  tests  in  the  houses  were  simulated  with  CONTAM93  by  including  a constant 
flow  element  in  the  door  of  each  house  and  adjusting  the  flow  until  a pressure  differences  of  4 
and  50  Pa  was  achieved.  The  airflow  rates  at  50  Pa  were  divided  by  the  interior  volumes  of  the 
houses  to  determine  the  50  Pa  air  change  rates,  and  the  4 Pa  flows  were  converted  to  effective 
leakage  areas  using  Equation  27  in  Chapter  23  of  ASHRAE  (3).  The  results  of  the  fan 
pressurization  simulations  are  shown  in  Table  1.  The  difference  between  the  Miami  and 
Minneapolis  houses  is  due  primarily  to  the  existence  of  the  basement  in  the  Minneapolis  houses. 
In  terms  of  both  measures  of  airtightness,  the  tight  houses  are  about  66%  tighter  than  the  houses 
of  typical  leakage. 


Table  1 - Fan  pressurization  simulation  results 


House 

achjo 

Leakage  area 

(hr-') 

(cm^) 

Typical  Miami  ranch 

13.2 

680 

Tight  Miami  ranch 

4.1 

220 

Typical  Minneapolis  ranch 

6.6 

720 

Tight  Minneapolis  ranch 

2.2 

230 

Typical  Miami  2 story 

12.9 

1,120 

Tight  Miami  2 story 

4.6 

390 

Typical  Minneapolis  2 story 

8.8 

1,170 

Tight  Minneapolis  2 story 

3.1 

410 

CONTAM93  was  used  to  calculate  whole  building  air  change  rates  for  wind  speeds  from  0 to  10 
m/s  and  indoor-outdoor  temperature  differences  from  -10  to  30  °C.  The  wind  direction  was  held 
constant  throughout  the  simulations.  These  simulations  were  performed  with  the  HVAC  systems 
both  on  and  off.  Whole  building  air  change  rates  were  calculated  by  adding  the  airflow  entering 
the  conditioned  space  of  the  house  through  all  leakage  paths.  The  results  of  these  airflow 
simulations  are  shown  in  Tables  2 through  9 for  the  system  off. 

Several  general  trends  are  shown  by  these  tables.  Using  'tight'  values  for  the  airflow  elements  vs. 
'typical'  or  best  estimate  values  reduced  the  whole  building  air  change  rate  by  up  to  a factor  of 
four  as  compared  to  a factor  of  three  for  the  fan  pressurization  results.  Also,  over  the  range 
considered  here,  the  wind  speed  had  a much  greater  impact  on  the  whole  building  air  change  rate 
than  the  temperature  difference.  However,  the  tight  airflow  elements  reduced  the  impact  of  the 
wind  speed  more  than  the  impact  of  the  temperature  difference. 
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Table  2 - Whole  house  airchange  rate  for typicalMiamiranch  house  (ach) 


Tin  - Tout  (K) 

-10 

-5 

0 

5 

10 

15 

20 

25 

30 

Wind  speed  (m/s) 
0 

0.33 

0.21 

0.00 

0.22 

0.35 

0.46 

0.57 

0.67 

0.76 

2 

0.40 

0.32 

0.33 

0.38 

0.47 

0.54 

0.65 

0.74 

0.84 

4 

0.75 

0.78 

0.82 

0.85 

0.89 

0.94 

1.00 

1.08 

1.15 

6 

1.31 

1.34 

1.38 

1.42 

1.46 

1.50 

1.54 

1.61 

1.67 

8 

1.92 

1.96 

2.01 

2.06 

2.11 

2.16 

2.21 

2.27 

2.33 

10 

2.57 

2.63 

2.69 

2.75 

2.81 

2.87 

2.94 

3.01 

3.08 

Table  3 - Whole  house  air  change  rate  for  tight  Miamiranch  house  (ach) 


Tin  - Tout  (K) 

-10 

-5 

0 

5 

10 

15 

20 

25 

30 

Wind  speed  (m/s) 
0 

0.10 

0.07 

0.00 

0.07 

0.11 

0.14 

0.17 

0.20 

0.23 

2 

0.11 

0.09 

0.08 

0.10 

0.14 

0.17 

0.20 

0.23 

0.26 

4 

0.18 

0.18 

0.19 

0.21 

0.22 

0.24 

0.26 

0.28 

0.31 

6 

0.30 

0.31 

0.32 

0.33 

0.34 

0.36 

0.38 

0.39 

0.42 

8 

0.44 

0.46 

0.47 

0.48 

0.49 

0.51 

0.53 

0.54 

0.57 

10 

0.60 

0.61 

0.63 

0.64 

0.65 

0.67 

0.69 

0.71 

0.73 

Table  4 - Whole  house  airchange  rate  for  typical  Minneapolis  ranch  house  (ach) 


Tin  - Tout  (K) 

-10 

-5 

0 

5 

10 

15 

20 

25 

30 

Wind  speed  (m/s) 
0 

0.25 

0.16 

0.00 

0.16 

0.26 

0.34 

0.42 

0.49 

0.56 

2 

0.29 

0.23 

0.18 

0.23 

0.31 

0.39 

0.46 

0.53 

0.59 

4 

0.45 

0.41 

0.44 

0.47 

0.50 

0.54 

0.59 

0.64 

0.69 

6 

0.69 

0.72 

0.75 

0.78 

0.81 

0.83 

0.87 

0.91 

0.95 

8 

1.03 

1.06 

1.09 

1.12 

1.16 

L19 

L22 

1.26 

1.29 

10 

1.39 

1.43 

1.46 

1.50 

1.54 

1.57 

1.62 

1.66 

1.70 

Table  5 - Whole  house  airchange  rate  for  tight  Minneapolis  ranch  house  (ach) 


Tin  - Tout  (K) 

-10 

-5 

0 

5 

10 

15 

20 

25 

30 

Wind  speed  (m/s) 
0 

0.09 

0.06 

0.00 

0.06 

0.10 

0.13 

0.16 

0.19 

0.21 

2 

0.09 

0.07 

0.04 

0.07 

0.11 

0.14 

0.17 

0.20 

0.22 

4 

0.13 

0.10 

0.11 

0.12 

0.14 

0.16 

0.19 

0.21 

0.24 

6 

0.17 

0.18 

0.19 

0.19 

0.21 

0.22 

0.23 

0.25 

0.27 

8 

0.25 

0.26 

0.27 

0.28 

0.29 

0.30 

0.31 

0.33 

0.34 

10 

0.34 

0.35 

0.36 

0.37 

0.38 

0.39 

0.41 

0.42 

0.44 
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Table  6 - Whole  house  air  change  rate  for  typical  Miami  2 story  house  (ach) 


Tin  - Tout  (K) 

-10 

-5 

0 

5 

10 

15 

20 

25 

30 

Wmd  speed  (m/s) 
0 

0.38 

0.24 

0.00 

0.25 

0.40 

0.53 

0.64 

0.76 

0.87 

2 

0.44 

0.34 

0.36 

0.42 

0.51 

0.62 

0.72 

0.81 

0.91 

4 

0.82 

0.86 

0.89 

0.93 

0.96 

1.02 

1.08 

1.15 

1.21 

6 

1.43 

1.47 

1.51 

1.55 

1.60 

1.64 

1.68 

1.74 

1.80 

8 

2.10 

2.15 

2.20 

2.25 

2.30 

2.36 

2.41 

2.47 

2.53 

10 

2.82 

2.88 

2.94 

3.01 

3.07 

3.14 

3.21 

3.28 

3.35 

Table  7 - Whole  house  air  change  rate  for  tight  Miami  2 story  house  (ach) 


Tin  - Tout  (K) 

-10 

-5 

0 

5 

10 

15 

20 

25 

30 

Wmd  speed  (m/s) 
0 

0.13 

0.08 

0.00 

0.09 

0.14 

0.18 

0.22 

0.26 

0.30 

2 

0.14 

0.10 

0.09 

0.12 

0.16 

0.21 

0.25 

0.28 

0.32 

4 

0.20 

0.22 

0.23 

0.24 

0.26 

0.28 

0.30 

0.34 

0.38 

6 

0.36 

0.37 

0.38 

0.40 

0.41 

0.43 

0.44 

0.47 

0.49 

8 

0.53 

0.54 

0.56 

0.57 

0.59 

0.60 

0.62 

0.64 

0.66 

10 

0.71 

0.73 

0.75 

0.76 

0.78 

0.80 

0.82 

0.84 

0.86 

Table  8 - Whole  house  air  change  rate  for  typicalMinneapolis  2 story  house  (ach) 


Tin  - Tout  (K) 

-10 

-5 

0 

5 

10 

15 

20 

25 

30 

Wmd  speed  (m/s) 
0 

0.25 

0.15 

0.00 

0.17 

0.27 

0.35 

0.43 

0.50 

0.58 

2 

0.30 

0.24 

0.25 

0.28 

0.34 

0.42 

0.48 

0.54 

0.61 

4 

0.57 

0.60 

0.62 

0.64 

0.66 

0.70 

0.74 

0.78 

0.83 

6 

0.99 

1.02 

1.05 

1.08 

1.10 

1.13 

1.16 

1.20 

1.24 

8 

1.46 

1.49 

1.52 

1.56 

1.60 

1.63 

1.67 

1.71 

1.75 

10 

1.95 

2.00 

2.04 

2.08 

2.12 

2.17 

2.22 

2.27 

2.32 

Table  9 - Whole  house  air  change  rate  fortight  Minneapolis  2 story  house  (ach) 


Tin  - Tout  (K) 

-10 

-5 

0 

5 

10 

15 

20 

25 

30 

Wmd  speed  (m/s) 
0 

0.09 

0.06 

0.00 

0.06 

0.09 

0.12 

0.15 

0.18 

0.20 

2 

0.10 

0.07 

0.06 

0.08 

0.11 

0.14 

0.17 

0.19 

0.21 

4 

0.14 

0.15 

0.16 

0.17 

0.18 

0.19 

0.21 

0.23 

0.26 

6 

0.25 

0.26 

0.27 

0.28 

0.29 

0.30 

0.31 

0.33 

0.34 

8 

0.37 

0.38 

0.39 

0.40 

0.41 

0.42 

0.44 

0.45 

0.47 

10 

0.50 

0.51 

0.53 

0.54 

0.55 

0.56 

0.58 

0.59 

0.61 
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Tables  10  through  17  present  the  results  of  the  airflow  simulations  with  the  HVAC  system  on. 
Operation  of  the  HVAC  system  increased  the  building  air  change  rate  as  much  as  0.3 1 ach  at  zero 
wind  speed  and  temperature  difference  due  to  supply  duct  leakage  in  the  attic.  The  effect  of  the 
system  fan  was  less  than  0.07  ach  at  high  wind  speeds  (>  4 m/s)  and  temperature  differences  (> 

10  °C). 


Table  10  - Wiole  house  air  change  rate  for  typical  Miami  ranch  house  with  system  on  (ach) 


Tm-Tout(R) 

-10 

0 

5 

10 

B 

20 

■ 25 

30 

Wind  speed  (m/s) 
0 

0.45 

0.38 

0.31 

0.39 

0.52 

0.63 

0.73 

0.83 

0.93 

2 

0.59 

0.52 

0.41 

0.50 

0.63 

0.74 

0.84 

0.93 

L03 

4 

0.86 

0.81 

0.85 

0.89 

0.95 

L02 

LIO 

L17 

L24 

6 

L34 

L37 

L41 

L45 

L49 

L55 

L61 

L67 

L73 

8 

L95 

L99 

2.04 

2.09 

2.14 

2.19 

225 

230 

2.38 

10 

2.60 

2.66 

2.72 

2.78 

2.84 

2.91 

2.97 

3.04 

3.11 

Table  11-  Wiole  house  air  change  rate  for  tight  Miami  ranch  house  with  system  on  (ach) 


tin  - lout  (K) 

-10 

-5 

0 

5 

10 

15 

20 

25 

30 

Wind  speed  (m/s) 
0 

029 

0.29 

0.30 

030 

0.30 

0.31 

031 

033 

037 

2 

0.30 

0.30 

0.30 

0.30 

030 

0.31 

032 

0.36 

039 

4 

0.37 

036 

0.34 

0.33 

031 

0.32 

0.36 

0.41 

0.44 

6 

0.46 

0.45 

0.44 

0.43 

0.41 

0.40 

0.42 

0.46 

0.49 

8 

0.56 

0.55 

0.53 

0.52 

0.53 

0.55 

0.57 

0.58 

0.61 

10 

0.65 

0.64 

0.64 

0.66 

0.69 

0.71 

0.73 

0.75 

0.77 

Table  12  - Wiole  house  air  change  rate  for  typical  Minneapolis  ranch  house  with  system  on  (ach) 


Tm-T6iit(Rr 

-10 

0 

5 

K5 

E 

20 

25 

30 

Wind  speed  (m/s) 
0 

0.24 

0.16 

0.00 

0.15 

0.25 

0.33 

0.41 

0.48 

0.55 

2 

0.28 

022 

0.18 

022 

0.30 

0.38 

0.45 

0.52 

0.58 

4 

0.44 

0.41 

0.44 

0.47 

0.50 

0.53 

0.59 

0.64 

0.68 

6 

0.69 

0.72 

0.75 

0.78 

0.80 

0.83 

0.87 

0.91 

0.94 

8 

1.03 

L06 

L09 

L12 

L16 

L19 

L22 

126 

129 

10 

L39 

L42 

L46 

L50 

L53 

L57 

L61 

L65 

L70 

Table  B - Wiole  house  air  change  rate  for  tight  Minneapolis  ranch  house  with  system  on  (ach) 


Im  - Tout  (K) 

-10 

— :5 — 

0 

— 5 — 

ID 

— E — 

■'^ir  ■ 

—1^ 

Wind  speed  (m/s) 
0 

0.09 

0.05 

0.00 

0.05 

0.09 

0.2 

0.15 

0.18 

021 

2 

0.08 

0.06 

0.04 

0.06 

0.10 

O.B 

0.16 

0.19 

021 

4 

0.12 

0.10 

0.11 

0.12 

O.B 

0.15 

0.18 

021 

023 

6 

0.17 

0.18 

0.18 

0.19 

020 

022 

023 

025 

026 

8 

025 

026 

027 

028 

029 

030 

0.31 

032 

034 

10 

0.34 

0.35 

036 

0.37 

038 

0.39 

0.40 

0.42 

0.43 
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Table  14  - Wiole  house  air  change  rate  for  typical  Miami  2 stoiy  house  with  system  on  (ach) 


Tm  - Tout  (K) 

-10 

0 

3 

ID — 

13 

20 

23 

30 

Wnd  speed  (m/s) 
0 

0.38 

024 

0.00 

025 

0.40 

033 

0.64 

0.76 

0.87 

2 

0.44 

034 

036 

0.42 

031 

0.62 

0.72 

0.81 

0.91 

4 

0.82 

0.86 

0.89 

0.93 

0.96 

L02 

L08 

L15 

121 

6 

L43 

L47 

152 

L56 

L60 

L64 

L68 

L74 

L80 

8 

2.10 

2.15 

220 

225 

231 

236 

2.41 

2.47 

2.53 

10 

2.82 

2.88 

2.94 

3.00 

3.07 

3.14 

321 

328 

335 

Table  15  - Wiole  house  air  change  rate  for  tight  Miami  2 story  house  with  system  on  (ach) 


Tm  - lout  (K) 

-10 

-5 

D 

3 

ID 

13 

20 

'■■23 

Wnd  speed  (m/s) 
0 

0.13 

0.08 

0.00 

0.09 

0.14 

0.18 

022 

026 

030 

2 

0.14 

0.10 

0.09 

0.2 

0.16 

021 

025 

028 

032 

4 

020 

022 

023 

024 

026 

027 

030 

034 

038 

6 

036 

037 

038 

039 

0.41 

0.43 

0.44 

0.47 

0.49 

8 

0.53 

0.54 

0.56 

0.57 

0.59 

0.60 

0.62 

0.64 

0.66 

10 

0.71 

0.73 

0.75 

0.76 

0.78 

0.80 

0.82 

0.84 

0.87 

Table  16  - \^hole  house  air  change  rate  for  typical  Minneapolis  2 story  house  with  system  on  (ach) 


rin-Tout(R) 

-10 

-5 

0 

5 

10 

13 

20 

23 

30 

Wnd  speed  (m/s) 
0 

025 

0.16 

0.01 

0.17 

027 

0.35 

0.43 

0.51 

0.58 

2 

031 

024 

025 

029 

035 

0.42 

0.48 

0.55 

0.61 

4 

0.57 

0.60 

0.62 

0.64 

0.66 

0.70 

0.74 

0.79 

0.83 

6 

0.99 

L02 

L05 

L08 

111 

L13 

L16 

120 

124 

8 

L46 

L49 

L53 

L56 

L60 

L63 

L67 

L71 

L75 

10 

L95 

2.00 

2.04 

2.08 

2.B 

2.17 

222 

227 

232 

Table  17  - Wide  house  air  change  rate  for  tight  Mirmeapolis  2 story  house  with  system  on  (ach) 


Tm  - Tout  (K) 

-10 

-5 

0 

5 

10 

15 

20 

23 

30 

Wnd  speed  (m/s) 
0 

0.09 

0.06 

0.00 

0.06 

0.09 

0.12 

0.15 

0.18 

020 

2 

0.10 

0.07 

0.06 

0.08 

0.11 

0.14 

0.17 

0.19 

021 

4 

0.15 

0.15 

0.16 

0.17 

0.18 

0.19 

021 

023 

026 

6 

026 

026 

027 

028 

029 

030 

031 

033 

034 

8 

038 

038 

039 

0.40 

0.41 

0.42 

0.44 

0.45 

0.47 

10 

0.50 

0.51 

0.53 

034 

0.55 

036 

0.58 

0.59 

0.61 
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Appendix  B Baseline  and  Preliminary  Simulation  Results 

Tables  la  through  24e  of  Appendix  B summarize  the  results  of  all  24  baseline  simulations. 

Tables  25a  through  27e  summarize  the  results  of  the  3 preliminary  simulations  of  the  lAQ 
control  retrofits.  Tables  la  through  27a  show  the  overall  peak  concentrations  (excluding  the 
basement,  attic,  garage  and  closet  zones),  the  location  of  that  overall  peak,  and  the  whole  house 
24-hour  average  concentrations  (excluding  the  basement,  garage,  and  attic  zones).  Tables  lb 
through  27b  show  the  individual  zone  peak  concentrations  for  the  main  living  space  zones. 

Tables  Ic  through  27c  show  the  individual  zone  24-hour  average  concentrations.  Tables  Id 
through  27d  show  the  individual  zone  4-hour  average  concentrations.  The  4-hour  average  was 
calculated  for  the  VOC  burst  sources  from  7 p.m.  to  1 1 p.m.,  for  the  oven  from  6 p.m.  to  10  p.m., 
and  for  the  heater  from  7 am  to  1 1 am.  No  4-hour  average  was  calculated  for  either  the  floor 
VOC  source  or  the  outdoor  air  pollutants.  Tables  le  through  27e  show  the  individual  zone 
1-hour  average  CO  concentrations.  The  1-hour  average  was  calculated  for  the  oven  from  7 p.m. 
to  8 p.m.  and  for  the  heater  from  9 am  to  10  am. 
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Table  la  - SIMIFLC  overall  24-hr  avg  concentrations 

VOCl  VOC2  VOC3  VOC4  VOC5  VOC6  VOC7  VOC8  VOC9  COH  N02.1  PART.l  C0.2  NOZ2  PART.2  COS  NOZ3  PART.3 

(Hg/m’)  (Mg/m’)  (Mg/m’)  (Mg/m’)  (Mg/m’)  (Mg/m’)  (Mg/m’)  (Mg/m’)  (Mg/m’)  ppm  ppm  (pg/m’)  ppm  ppm  (pg/m’)  (PPM)  (PPM)  (pg/m’) 

24  hr  avg  265  6067  218  183  117  225  219  217  205  2.75  0.026  10,86  1.60  0.008  10,75  6.77  0,079  14,63 


Table  2a  - SIMIFLM  overall  24-hr  avg  concentrations 

VOCl  VOC2  VOC3  VOC4  VOC5  VOC6  VOC7  VOC8  VCM:9  CO.l  N02.1  PART.l  C0.3  N02.3  PART.3 

(Mg/m»)  (Mg/m»)  (Mg/m’)  (^g/m’)  (Mg/m’)  (Mg/m’)  (Mg/m’)  (Mg/m’)  (gg/m’)  ppm  ppm  (pg/m’)  ppm  ppm  (pg/m’) 

24hravg  235  4899  165  152  98  206  189  218  183  2.50  0,025  11,41  6.68  0,086  18.45 


Table  3a  - SIMIFLH  overall  24-hr  avg  concentrations 

VOCl  VOC2  VOC3  VOC4  VOC5  VOC6  VOC7  VOC8  VOC9  CO.l  N02.1  PART.l  C0.3  N02.3  PART.3 

(Hg/m’)  (Mg/m’)  (Mg/m’)  (Mg/m’)  (Mg/m’)  (Mg/m’)  (Mg/m’)  (Mg/m’)  (Mg/m’)  ppm  ppm  (pg/m’)  ppm  ppm  (pg/m’) 

24hravg  319  6289  195  199  191  196  201  198  198  2,54  0,023  8.93  6.97  0.078  7,52 


Table  4a  - SIM  1 FTC  overall  24-hr  avg  concentrations 

VOCl  VOC2  VOC3  VOC4  VOC5  VOC6  VOC7  VOC8  VOC9  CO.l  N02.1  PART.l  C0.2  N02.2  PART.2  C0.3  N02.3  PART.3 

(Hg/m>)  (Mg/m»)  (ng/m^)  (pg/m^)  (ng/m^)  (pg/m*)  (jig/m^)  (ng/m^)  (ng/m’)  ppm  ppm  (pg/m’)  ppm  ppm  (pg/m^)  (PPM)  (PPM)  (pg/m’) 

24hravg  768  19246  451  402  181  460  472  453  449  4.63  0.025  7.95  1.66  0.003  7.80  6.66  0.032  5.02 


Table  5a  - SIM  IFTM  overall  24-hr  avg  concentrations 

VOCI  VOC2  VOC3  VOC4  VOC5  VOC6  VOC7  VOC8  VOC9  CO.l  N02.1  PART.l  C0.3  N02.3  PART.3 

(Mg/m»)  (fjg/m*)  (Mg/m*)  (Mg/m*)  (Mg/m*)  (Mg/m*)  (Mg/m*)  (Mg/m*)  (Mg/m*)  ppm  ppm  (pg/m*)  ppm  ppm  (pg/m*) 

24hravg  682  21165  485  407  136  525  507  505  510  4.18  0.025  9.23  6,94  0,039  7.32 


Table  6a  - SIMIFTH  overall  24-hr  avg  concentrations 

VOCl  VOC2  VOC3  VOC4  VOC5  VOC6  VOC7  VOC8  VOC9  CO.l  N02.1  PART.l  C0.3  N02.3  PART.3 

(Mg/m^)  (Mg/m’)  (pg/m’)  (pg/m’)  (Mg/m’)  (Mg/m’)  (gg/m’)  (Mg/m’)  (gg/m’)  ppm  ppm  (pg/m’)  ppm  ppm  (pg/m’) 

24hravg  465  9517  239  243  507  237  249  237  240  2,91  0.022  7.55  6.99  0,056  4,57 


Table  7a  - SIM  1 MLC  overall  24-hr  avR  concentrations 

VOCl  VOC2  VOC3  VOC4  VOC5  VOC6  VOC7  VOC8  VOC9  CO.l  N02.1  PART.l  C0.2  N02.2  PART.2  C0.3  N02.3  PART.3 

(Hg/m»)  (Mg/m*)  (Mg/m*)  (Mg/m*)  (Mg/m*)  (Mg/m*)  (Mg/m*)  (Mg/m*)  (Mg/m*)  (PPM)  (PPM)  (Mg/m*)  (PPM)  (PPM)  (Mg/m*)  (PPM)  (PPM)  (Mg/m*) 

24hravg  182  2767  144  140  127  145  143  364  144  2,01 0,020  10,23  1,84  0.018  10.58  6.70  0,102  14,46 


Table  8a  - SIMIMLM  overall  24-hr  avg  concentrations 

VOCI  VOC2  VOC3  VCK:4  VOC5  VOC6  VOC7  VOC8  VOC9  CO.l  N02.1  PART.l  C0.2  N02.2  PART.2  C0.3  N02.3  PART.3 

(Hg/m^)  (Mg/m^)  (ng/m*)  (pg/m*)  (Mg/m*)  (Mg/m*)  (Mg/m*)  (Mg/m*)  (Mg/m*)  (PPM)  (PPM)  (Mg/m*)  (PPM)  (PPM)  (Mg/m*)  (PPM)  (PPM)  (Mg/m*) 

24hravg  265  3489  133  146  176  145  133  325  153  1.99  0.020  10,85  1.93  0.018  11.32  6.69  0.096  18.68 


Table  9a  - SIMIMLH  overall  24-hr  avg  concentrations 

VOCl  VOC2  VOC3  VOC4  VOC5  VOC6  VOC7  VOC8  VOC9  CO.l  N02.1  PART.  I CO.  3 N02.3  PART.3 

(fig/m^)  (Mg/m’)  (Mg/m’)  (Mg/m’)  (Mg/m’)  (Mg/m’)  (Mg/m’)  (Mg/m’)  (Mg/m’)  (PPM)  (PPM)  (Mg/m’)  (PPM)  (PPM)  (Mg/m’) 

24  hr  avg  243  3171  132  137  131  136  128  333  141  1,94  0,020  10,59  6.67  0.099  23.44 


Table  10a  - SIMIMTC  overall  24-hr  avg  concentrations 

VOCl  VOC2  VOC3  VOC4  VOC5  VOC6  VOC7  VOC8  VOC9  CO.l  N02.1  PART.l  C0.2  N02.2  PART.2  C0.3  N02.3  PART.3 

(Mg/m*)  (ng/m^)  (Mg/m’)  (ng/mS)  (lig/m^)  (ng/m^)  (^g/m’)  (jig/m»)  (ng/m*)  (PPM)  (PPM)  (pg/m^)  (PPM)  (PPM)  (pg/m^)  (PPM)  (PPM)  (pg/m^) 

24hravg  213  5161  185  186  266  187  183  728  188  2,79  0.018  6.83  2.43  0.016  7.62  6.66  0.045  4.96 


Table  1 la  - SIM  I MTM  overall  24-hr  avg  concentrations 

VOCl  VOC2  VOC3  VOC4  VOC5  VOC6  VOC7  VOC8  VOC9  CO.  I N02.1  PART.  I C0.2  N02.2  PART.2  C0.3  N02.3  PART.3 

(Hg/m^)  (Mg/in’)  (jig/m’)  (ng/m^)  (jig/in’)  (ng/m’)  (ng/m>)  (^lg/m^)  (ng/nP)  (PPM)  (PPM)  (ng/nP)  (PPM)  (PPM)  (Mg/m’)  (PPM)  (PPM)  (ng/m’) 

24hravg  526  10674  235  257  295  254  231  1144  270  2.85  0.021  7.84  2.79  0.012  8.91  6.72  0.041  6.67 
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Table  12a  - SIMIMTH  overall  24-hr  avg  concentrations 

VOCl  VOC2  VOC3  VOC4  VOC5  VOC6  VOC7  VOC8  VOC9  CO.l  N02.1  PART.l  C0.3  N02.3  PART.3 

(pg/m^)  (pg/m^)  (pg/m^)  (pg/m^)  (pg/m*)  (pg/m^)  (pg/m’)  (pg/m*)  (pg/m^)  (PPM)  (PPM)  (pg/m^)  (PPM)  (PPM)  (pg/m^) 

24  hr  avg  446  9448  227  232  205  232  223  1025  238  2.80  0,019  7.38  6.88  0.044  9.09 


Table  13a  - SIM2FLC  overall  24-hr  avg  concentrations 

VOCl  VOC2  VOC3  VOC4  VOC5  VOC6  VOC7  VOC8  VOC9  CO.l  N02.1  PART.!  C0.2  N02.2  PART.2  C0.3  N02.3  PART.3 

(pg/m>)  (Mg/nP)  (Mg/m’)  (Mg/m*)  (Mg/m’)  (Mg/m’)  (Mg/m’)  (Mg/nP)  (Mg/m’)  (PPM)  (PPM)  (Mg/m’)  (PPM)  (PPM)  (Mg/m’)  (PPM)  (PPM)  (Mg/m’) 

24hfavg  813  9167  127  165  149  178  144  121  157  2.26  0.019  10.87  I.S9  0,009  10,81  6.68  0,086  15.79 


Table  14a  - S1M2FLM  overall  24-hr  avg  concentrations 

VOCI  VOC2  VOC3  VOC4  VOC5  VOC6  VOC7  VOC8  VOC9  CO.l  N02.1  PARTI  C0.3  N02.3  PART.3 

((jg/nP)  (jig/nP)  ((ig/nP)  (ng/nP)  (pg/nP)  (jig/nP)  ((ig/nP)  (pg/nP)  (Mg/nP)  (PPM)  (PPM)  (pg/m’)  (PPM)  (PPM)  (pg/nP) 

24  hr  avg  681  7974  135  137  134  I4S  167  99  IS3  1.94  0 017  11.22  6.60  0.092  20.38 
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Table  14e  - SIM2FLM  zone  1-hr  avg  concentrations  LEGEND 


Table  15a  - SIM2FLH  overall  24-hr  avg  concentrations 

VOCI  VOC2  VOC3  VOC4  VOC5  VOC6  VOC7  VOC8  VOC9  CO.l  N02.1  PART.I  C0.3  N()2.3  PARr.3 

((ig/nP)  (jig/nP)  ((ig/m’)  (Mg/ni*)  (pg/nP)  (Mg/in’)  (Mg/m*)  (pg/nP)  (pg/nP)  (PPM)  (PPM)  (pg/m>)  (PPM)  (Pl’M)  (pg/nP) 

24  hr  avg  1441  13042  186  178  178  181  189  102  185  2,37  0.016  8,49  6.93  0.063  7.92 


Table  16b  - SIM2FTC  zone  peak  concentrations 
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Table  16e  - SIM2FTC  zone  l-hr  avg  concentrations  LEGEND 


Table  17a  - S1M2FTM  overall  24-hr  avg  conccnlrations 

VOCI  VOC2  VOC3  VOC4  VOC5  VOC6  VOC7  VOC8  VOC9  CO.l  N02.1  PARTI  C0.3  N02.3  PART.3 

(pg/m»)  (Mg/nP)  (Mg/nP)  (ng/nP)  (ng/nP)  (pg/nP)  (pg/nP)  (pg/nP)  (pg/nP)  (PPM)  (PPM)  (pg/m>)  (PPM)  (PPM)  (pg/nP) 

24  Ilf  avg  1828  23498  213  222  228  211  282  106  272 2.S9  0.015  9.13  6.78  0,044  8.64 
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Table  1 8a  - S1M2FTH  overall  24-hr  avg  concentrations 

VOCl  VOC2  VOC3  VOC4  VOC5  VOC6  VOC7  VOC8  VOC9  CO.l  N02.1  PART.l  C0.3  N02.3  PART.3 

(Hg/nP)  ((ig/nP)  (ng/nP)  (pg/nP)  (pg/nP)  (pg/nP)  (pg/irP)  (pg/itP)  (pg/itP)  (PPM)  (PPM)  (pg/m’)  (PPM)  (PPM)  (pg/m’) 

24hfavB  3979  39385  370  359  360  359  369  122  367  3.92  0.015  4,94  7.00  0.024  2,89 


Tabic  19a  - SIM2MLC  overall  24-hr  avg  concentrations 

VOCI  VOC2  VOC3  VOC4  VOC5  VCX:6  VOC7  VOC8  VOC9  CO.  I N02.1  PART.l  C0.2  N02.2  PART.2  C0.3  N02.3  PART.3 

(Hg/m*)  ((ig/m*)  (Mg/m’)  (pg/m’)  (pg/™’)  (pg/nP)  (pg/m’)  (pg/m>)  (pg/m>)  (PPM)  (PPM)  (pg/nP)  (PPM)  (PPM)  (pg/iiP)  (I’PM)  (PPM)  (pg/m>) 

24lifavg  248  3586  115  126  124  128  116  141  121  1.86  0.018  10.67  1.76  0.019  10,92  6,73  0,114  16.48 


Table  20a  - S1M2MLM  overall  24-hr  avg  concentrations 

VOCl  VOC2  VOC3  VOC4  VOC5  VOC6  VOC7  VOC8  VOC9  CO.I  N02.1  PARTI  C0.2  N02.2  PART.2  C0.3  N02.3  PART.3 

(Hg/m>)  (Mg/nP)  (Mg/m>)  (jig/np)  (|ig/ni>)  (ng/m>)  ((ig/m’)  (ng/m’)  (ng/m’)  (PPM)  (PPM)  (pg/m’)  (PPM)  (PPM)  ((ig/m*)  (PPM)  (PPM)  (ng/m>) 

24tuavR  269  6490  129  125  116  136  127  159  126  1,88  0.017  10.76  1.82  0,017  H.I2  6.60  0.095  19.19 
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Table  2 la  - S1M2MLH  overall  24-hr  avg  concenlralions 

VOCI  VOC2  VOC3  vcx:4  VOC5  voce  VOC7  VOC8  VOC9  CO.l  N02.1  PARTI  C0.3  N02.3  PART.3 

(pg/nP)  (pg/nP)  (pg/ntP)  (pg/nP)  (pg/nP)  (pg/nP)  (pg/nP)  (pg/nP)  (pg/nP)  (PPM)  (PPM)  (pg/nP)  (PPM)  (PPM)  (pg/nP) 

24hfavg  284  7364  123  I2S  134  127  128  132  123  1.87  0.017  10.73  6.71  0.097  17.87 


Table  22b  - S1M2MTC  zone  peak  concentrations 
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Tabic  23a  - SIM2MTM  overall  24-hr  avg  concenlralions 

VOCI  VOC2  VOC3  VOC4  VOC5  VOC6  VOC7  VOC8  VOC9  CO.  I N02.1  PART.l  C0.2  N02.2  PART.2  C0.3  N02.3  PART.3 

(Mg/m>)  (ng/m>)  ((ig/m’)  (pg/m>)  (ng/m’)  ((ig/ni»)  (jig/m’)  (ng/m’)  ((ig/nP)  (PPM)  (PPM)  ((ig/m’)  (PPM)  (PPM)  ((ig/nP)  (PPM)  (PPM)  (Mg/m>) 

24hravg  862  23366  184  202  220  216  207  261  195  2.55  0.016  7.96  2,44  0,012  8.85  6.58  0.040  6.77 


Tabic  24a  - S1M2MTH  overall  24-hr  avg  concentrations 

VOCI  VOC2  VOC3  VOC4  VOC5  VOC6  VOC7  VOC8  VOC9  CO.l  N02.1  PART.I  CO,3  N02.3  PART.3 

(pg/m’)  (pg/nP)  {pg/m>)  (pg/nP)  (pg/irp)  (pg/m’)  (pg/m’)  (pg/m>)  (pg/m’)  (PPM)  (PPM)  (pg/m’)  (PPM)  (PPM)  (pg/m>) 

24hravg  802  21581  184  205  218  210  200  232  197  2.60  O.OIS  7.78  6.83  0.041  6.13 
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Table  24e  - S1M2MTH  zone  1-hr  avg  concentralions  LEGEND 


Table  25a  - SIM  IFLCF  overall  24-hr  avg  concentrations 

VOCl  VOC2  VOC3  VOC4  VOC5  VOC6  VOC7  VOC8  VOC9  CO.l  N02.1  PART.l  C0.2  N02.2  PART.2  C0.3  N02.3  PART.3 

(Hg/m’)  (Mg/m^)  (Mg/m^)  (Mg/m^)  (ng/m*)  (Mg/m’)  (pg/m’)  (pg/m’)  (pg/m’)  ppm  ppm  (pg/m’)  ppm  ppm  (pg/m’)  (PPM)  (PPM)  (pg/m’) 

24hravg  265  6067  218  183  117  225  219  217  205  2.75  0.026  9.26  1.60  0.008  9.17  6.77  0,079  14.52 
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Table  27a  - SIMIFLCHO  overall  24-hr  avg  concentrations 

VOCl  VOC2  VOC3  VOC4  VOC5  VOC6  VOC7  VOC8  VOC9  CO.l  N02.1  PART.l  C0.2  N02.2  PART.2  C0.3  N02.3  PART.3 

(Hg/m»)  (Mg/m^)  ((ig/m*)  (Mg/m’)  (Mg/m’)  (Mg/m’)  (Mg/m’)  (Mg/m’)  (Mg/m’)  ppm  ppm  (pg/m’)  ppm  ppm  (pg/m’)  (PPM)  (PPM)  (pg/m’) 

24hravg  253  5717  212  178  112  220  214  211  200  2.71  0.026  10,97  1.60  0.008  10.87  6.78  0,082  14.38 


